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1.0 


INTRODUCTION 


1.1  Background 

Extended,  randomly  time  varying  radar  targets,  clutter,  and 
communication  channels  have  been  considered  in  the  past  [1-4].  An 
important  additional  variable  for  such  problems  is  considered  here,  namely, 
random  polarization  modulation  [5].  Sea  echo,  for  example,  exhibits 
random,  Doppler* dependent  polarization  effects  [6]. 

The  design  of  radar  or  communications  systems  for  randomly  time 
varying  targets  or  channels  has  typically  been  accomplished  by  using  one 
of  two  approaches:  (i)  a  likelihood  ratio  test  [4]  or  (ii)  implementation 
of  a  receiver  with  maximum  signal- to-interference  ratio  (SIR)  [7-11].  Both 
of  these  approaches  often  use  a  scattering  function  description  of  the 
target,  clutter,  and  channel. 

SIR  maximization  is  typically  more  straightforward,  does  not  depend 
on  a  Guassianity  assumption,  and  gives  usable  results.  Maximization  of 
SIR,  however,  does  not  necessarily  result  in  an  optimum  receiver,  i.e.» 
an  implementation  of  a  likelihood  ratio  test.  For  the  special  case  of  a 
known  signal  in  colored  noise,  both  SIR  maximization  and  a  likelihood 
ratio  test  yield  the  same  result  (a  whiten  and  match  filter),  but  this 
correspondence  does  not  always  hold. 

For  Gaussian  data,  a  likelihood  ratio  test  for  random,  extended 
targets  can  be  implemented  if  the  scattering  functions  of  target  and  clutter 
are  known,  since  the  covariance  functions  of  the  echoes  are  then  also 
known  [12,13].  There  are,  however,  no  straightforward  signal  design 
techniques  associated  with  the  likelihood  ratio  test,  and  it  is  difficult 
to  define  such  a  test  in  the  polarimetric  case  [5]. 
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The  SIR  maximization  technique  would  be  especially  attractive  if: 
(i)  it  could  be  modified  in  such  a  way  as  to  implement  a  Bayes  optimum 
(likelihood  ratio)  receiver,  (ii)  the  performance  of  such  a  receiver  could 
be  predicted,  and  (iii)  SIR  maximization  could  be  applied  to  the  design 
of  polarization -sensitive  radar  and  communication  systems.  This  report 
shows  that  all  three  goals  are  attainable,  and  gives  specific  examples. 

The  basic  problem  is  illustrated  in  Figure  1-1,  along  with  some 
of  the  notation  used  in  the  sequel.  The  problem  is  to  obtain  a  suitable 
representation  of  doubly  spread  target  and  clutter  (a  polarimetric 
scattering  function) ,  and  to  use  this  representation  to  obtain  optimum 
vertically  and  horizontally  polarized  signal  and  filter  functions  u1(t), 
Ujft),  f  1< t)  and  f2(t).  It  is  also  desirable  to  generalize  the  receiver 
in  Figure  1-1  to  implement  a  likelihood  ratio  test  and  to  predict  the 
performance  of  such  a  test. 

1.2  Overview 


Section  2  analyzes  the  relation  between  maximization  of  signal-to- 
interference  ratio  and  a  likelihood  ratio  test  for  discrimination  of  zero  mean 
random  signals  in  two  channels.  Section  3  establishes  the  dependence  of 
SIR  on  polarimetric  scattering  functions.  The  connection  between 
scattering  functions  and  a  tapped  delay  line  scattering  model  with 
randomly  time  varying  tap  weights  is  also  given  in  Section  3,  The  tap 
weights  describe  energy  coupling  from  one  polarization  channel  to  another. 
Section  4  introduces  physical  insight  and  specific  mathematical  models  into 
the  polarimetric  scattering  function  formulation  by  considering  planar 
point  targets  and  randomly  tilted  dipoles  with  a  restricted  maximum  tilt 
relative  to  vertical.  Section  5  gives  some  mathematical  details  necessary 
for  implementation  of  a  computer  program  for  polarimetric  SIR  optimization. 
Sections  5.3  and  5.4  can  be  skipped  if  the  reader  is  not  interested  in 
such  details.  Sections  6  and  7  describe  the  application  of  the  computer 


2 


program  to  some  specific  examples.  Analysis  of  the  SIR  expression  for 
detecting  a  distribution  of  planar  reflectors  in  a  background  of  randomly 
oriented  dipoles  (distributed  planar  target  in  chaff)  yields  a  new 
polarimetric  chaff  cancellation  method.  This  method  is  obtained  in 
Section  7.6.  A  disparity  between  the  SIR  maximization  criterion  and  the 
power  of  a  likelihood  ratio  test  appears  in  the  context  of  a  receiver  with 
multiple  orthogonal  filters,  which  converts  the  maximum  SIR  processor  to 
a  Bayes  optimum  processor.  This  disparity  between  SIR  and  power 
measures  is  analyzed  in  Section  7.7.  A  review  of  the  results  is  given 
in  Section  8.0. 
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Figure  l-l.  Polarimetric  Radar  Configuration 


2.0  SIR  MAXIMIZATION  AND  OPTIMUM  HYPOTHESIS  TESTS 


The  following  analysis  demonstrates  that  orthogonal  filters  which 
maximize  SIR  can  also  be  used  for  simultaneous  diagonalization  of  both 
signal  and  interference  covariance  matrices.  The  filter  responses  are 
independent  random  variables  under  either  hypothesis,  and  the  log 
likelihood  statistic  is  easily  obtained.  Receiver  performance  is  also  easily 
obtained  in  closed  form.  These  results  are  important  because  SIR 
maximization  yields  a  "best”  signal  as  well  as  an  appropriate  set  of 
orthogonal  filters.  A  direct  likelihood  ratio  approach  is  based  on  echo 
covariance  matrices  which  can  only  be  evaluated  after  the  signal  has 
been  specified. 


2, 1  Signal-to- Interference  Ratio  (SIR)  Maximization  and 
an  Eigenfunctioh  Equation 


The  signal  to  interference  ratio  is 


CO 

1 

1  f  rMx)  f(x)  dx|2 

signal 

i 

J  - 

SIR 


y*  r*(*i  nx) 


dx 


inter  ferance 


(2-1) 


where  r*(x)  is  the  conjugate- transpose  of  the  2  *  1  data  column  vector. 
This  data  vector  i#  composed  of  a  vertically  polarized  component,  r^(x), 
and  a  horizontally  polarized  component,  rgfx).  Thus 
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E 


r*fx)  «  [  r  *  (  x)  r  *  ( x)  ] 

where  r*(x)  is  the  conjugate  of  the  complex  scalar  time  function  r^(x). 
Similarly,  f(x)  represents  a  2  *  1  column  vector  representing  the  time 
function  f^(x)  that  is  to  be  correlated  with  the  vertically  polarized  data, 
and  another  function  f g(x)  that  is  to  be  correlated  with  horizontally 
polarized  data.  These  functions  are  illustrated  in  Figure  1-1. 


Using  the  identity 


00  ' 

i  2 

r*(x)  f  ( x)  dx| 

-00  i 


00  00 

f*(x)  E  {  r  <  x )  r*  <y) }  fCy)  dxdy 

•  00  .00 

00  00 

-//  f*(x)  C  (x,y)  f (y)  dxdy 

—  GJ  •CD 


we  have 


09  at 


//  £* ( x)  Cs  (x,y)  £{y)  dydx 


SIR 


-00 


If  f*(x)  C£  tx,y)  f C y )  dydx 


(2-2) 


(2-3) 


where  Cg(x.y)  is  the  signal  covariance  matrix.  The  interference 
covariance  matrix  is 

( x, y)  -  Cc(x,y)  *  (Nq/2)  I  6<x  -  y) 
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where  C^,  is  the  2x2  clutter  covariance  matrix  and  N  q  /  2  is  the  noise 
power  spectral  density. 

To  maximize  SIR,  one  can  maximize  the  numerator  of  (2-3)  with  a 
constraint  on  the  denominator.  If  the  denominator  is  constrained  to  equal 
one,  an  equivalent  problem  to  maximizing  (2-3)  is  to  maximize  the  functional 


Q( 


f)  =  =S  £(y)  dydx 


-  X, 


f* 


f*(x)  C  (x,y)  f ( y )  dyds 


ffc-'U)  Cg(x,y)  -  XQ  Cj  ( x ,  y ) 1  f(y)  dyd 


(2-5) 


where  unlabelled  integration  limits  are  (-»,  «)  and  XQ  is  any  positive 
constant.  Another  constraint  is  to  make  the  filter  energy  equal  to 
unity,  i.e. , 


Egy(f)  =  / f* (x)  f (x)  dx  n  1 


(2-6) 


The  energy  constraint  is  not  included  explicitly  in  (2-3)  because  SIR  is 
invariant  when  f(x)  is  multiplied  by  a  nonzero  scalar  constant.  Such 
a  constraint,  however,  can  easily  be  incorporated  in  a  computer  optimization 
algorithm.  If  Egy(f)  is  specified  as  in  (2-6),  then  the  problem  is  to 
maximize 


fi 


f * ( x)  g ( x)  dx 


(2-7) 
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when  (2-6)  holds  and  when 


£(x)  -y*[Ss(x,y)  -  xD£x(x/y) 


f  (y)  dy 


By  the  Schwarz  inequality,  (2-7)  is  maximized  when  f(x)  is 
square-integrable  as  in  (2-6)  and  when 

£(x)  =»  k  f  (X) 

i.e.,  when  g(x)  is  proportional  to  f(x) .  Substituting  (2-8)  into  (2-9) 
we  have 

y^iglx/y)  -  Cj (x,y)  f  (y)  dy  *  k  f(x) 

Equation  (2-10)  will  be  satisfied  if  f(y)  is  an  eigenfunction  of  both 
Cg(x.y)  and  C^x.y),  i.e.,  if 

f  £g(x,y)  f (y)  dy  »  Xg  f (x) 
f  Cj(x,y)  (y )  dy  ■  Xj  £  (x) 

A  more  general  solution  of  (2-10)  can  be  obtained  if 

f  £s(x,y)  jlw  dy "  + 

f  Cj. (x,y)  f(y>  dy  -  Xj  fix)  +  d(x) 


and 


and 


(2-8) 


(2-9) 


(2-10) 


(2-11) 


(2-12) 


(2-13) 


(2-14) 
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where  b(x)  and  d(x)  are  such  that  b(x)  -  X^  d(x)  is  zero  for  any  positive 
constant,  But  b(x)  equals  d(x)  for  any  XQ  only  if  b(x)  =  d(x)  =  0, 
and  we  are  back  to  (2-11)  and  (2-12), 


Combining  (2-11)  and  (2-12)  yields 


fcsU,y)  f(y)  dy  =  <XS/V  f=  l<y>  d V  •  (2*15) 

To  obtain  further  insight  into  (2-15),  we  can  define  an  Inverse  kernel  as 
follows.  If  Cj^fx.y)  is  the  inverse  kernel  of  Cj(x,y)  then 


/ 


c/u,*)  C_(z,y)  da  ®  I  6{x  -  y) 

o  I  <=I  a*  •* 


(2-18) 


Applying  this  definition  to  (2-15),  we  have 

csc2,y>  !<y>  <*y<**  -  (xs/Xj>  «<x>  .  (2-n) 

The  solution  f(x)  to  the  SIR  maximization  problem  is  then  an  eigenfunction  of 


C  ( x, y )  *» 


z)  Cg (z, y)  dz 


(2-18) 


with  eigenvalue  Xg/Xj.  The  significance  of  this  eigenvalue  emerges  when 
(2-11)  and  (2-12)  are  substituted  into  (2-3): 


sir  - 


Xs/r<x>  fix)  dx  ^ 

s 

TT 

X.  /£Mx)  f(x)  dx  I 


</f* 


(2-19) 


The  energy -constrained  filter  that  maximizes  SIR  is  thus  the  eigen¬ 
function  of  £(x,y)  in  (2-18)  with  largest  possible  eigenvalue.  This 
relation  does  not  imply  that  one  can  easily  obtain  the  best  f(x)  from  an 
eigenfunction  equation,  however,  since  C(x,y)  is  undefined  until  the 
best  signal,  u(x),  is  specified.  The  SIR  maximization  approach  allows  one 
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to  obtain  both  u>.x)  and  f(x).  To  obtain  the  best  signal  with  an  eigen¬ 
function  formulation,  different  signals  must  be  tried,  the  corresponding 
covariance  matrices  C(x,y)  must  be  calculated,  and  the  covariance  matrix 
with  largest  principal  eigenvalue  must  be  identified. 

If  the  unit  energy  filter  function  f  t(x)  that  maximizes  SIR  is  the 
piincipai  eigenfunction  of  C(x,y)  in  (2-18),  then  can  the  other  eigenfunctions 
of  C(x,y)  also  be  obtained  by  SIR  maximization?  Consider  another  filter 
function,  f2(x),  with  the  following  properties: 

(i)  f2(x)  is  orthogonal  to  f^x),  i.e., 

f £i(x)  £2{x>  dx  *  0  '  (2-20) 

(ii)  £2fx)  has  un*t  ener*y  as  ^  (2-6),  and 
(iii>  f2(x)  maximizes  SIR. 

From  the  analysis  in  (2- 5) -(2- la) ,  £2(x>  U»  an  eigenfunction  of  C(x,y)  in 
(2-18),  if  such  an  eigenfunction  satisfies  (2-20).  In  fact,  all  the  eigenvectors 
of  a  covariance  matrix  with  distinct  eigenvalues  are  orthogonal  ( 14] ,  and 
(2-20)  is  satisfied. 

All  the  eigenvectors  of  C(x,y)  can  thus  be  obtained  by  SIR 

maximization,  provided  that  each  new  filter  function  is  constrained  to  be 

orthogonal  to  those  found  previously.  If  the  maximum  possible  SIR  is 

obtained  for  each  filter  function,  then  (2-17)  implies  that  the  n  computed 

tH 

eigenfunction  will  have  the  n  largest  eigenvalue,  as  in  principal  component 
analysis  ( 15] . 

Maximization  of  SIR  under  ortho jonality  constraints  as  in  (2-20)  is 
particularly  straightforward  if  the  simplex  method  (16]  is  used.  If  the 
initial  simplex  is  constrained  to  a  subspace  of  R  4,  then  the  solution  will 
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be  constrained  to  the  same  subspace.  The  starting  points  for  SIR 
maximization  will  correspond  to  functions  that  are  orthogonal  to  previously 
determined  filters  if  Gram-Schmidt  orthogonalization  is  used.  After 
constraining  the  starting  points  (i.e.,  the  vertices  of  the  initial  simplex) 
to  be  orthogonal  to  previously  determined  filters,  the  simplex  algorithm 
can  be  run  without  further  modification. 


2.2  Significance  of  the  Relation  between  SIR  Maximization  and 
Eigenfunction  Analysis 

By  maximizing  signal-to-interference  ratio,  one  can  obtain  the 
eigenfunctions  and  eigenvalues  of  C(x,y)  in  (2-18).  The  eigenfunctions 
are  found  in  two  contexts  in  the  literature.  First,  they  are  the  best  set 
of  linear  discriminants  for  discriminating  between  two  zero-mean  Gaussian 
processes  [17,18,19].  Second,  they  can  be  used  for  simultaneous 
diagonalization  of  both  signal  and  interference  covariance  matrices  [  20] . 

If  the  data  are  projected  along  the  N  eigenvectors  of  C(x,y)  in  (2-18), 
the  resulting  projections  will  be  uncorrelated,  with  variances  U 

m  on  n**i 

for  the  signal  process  and  for  the  interference  process. 

This  observation  follows  from  the  fact  that  each  eigenfunction  must 
satisfy  both  (2-11)  and  (2-12).  Another  proof  is  given  in  [20]. 


Simultaneous  diagonalization  yields  a  likelihood  ratio  formulation 

involving  simple  operations  on  the  outputs  of  the  filters  that  maximize  SIR. 

th 

Under  hypothesis  the  output  of  the  n  niter  is 


(2-21) 


the  sum  of  the  variances  of  the  uncorrelated  echo  and  interference  processes. 
Under  HQ,  only  the  interference  is  present  and 

•  <2-“> 


ii 


Because  of  simultaneous  diagonalization  and  Gaussianity,  the  filter  outputs 
are  independent  under  both  (signal  +  interference)  and  Hq  (interference 
alone). 

The  likelihood  ratio  is  then 

N 

A(r)  *  II  A(r  )  (2-23) 

,  n 


A  A 


where  r  is  the  vector  of  N  filter  responses,  r^,  r2,  ...  ,  rN«  Since  rn 
is  complex,  we  have  [21,22,23] 


A(r  ) 
n 


t*,XSn  *  S.”’1  expI*i;n|J/asn  *  Xln)] 


(2-24) 


til 

The  log-likelihood  ratio  for  the  n  filter  output  is 


Urn)  -  -mu  ♦  USn/AInM  * 


Xsn/Xm 

LXSn  *  XIn 


r 

n ‘ 


(2-25) 


The  log-likelihood  ratio  can  be  written  strictly  in  terms  of  SIR 
if  the  filter  outputs  are  whitened  before  being  passed  through  a  likelihood 
ratio  test.  To  whiten  the  interference,  the  n  filter  output  is  multiplied 
by  Xm 


In  ’ 


r 

n 


-1/ 
I 


r 


n 


(2-26) 
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For  the  whitened  filter  outputs,  we  have 


(  2  ) 

E  |  | r n l  i  H0  =  1  ,  n  *  l,  2,  ....  R 


(2-27) 


and 


E  Jlcnl2  l*x|  '  (Xsn/XIn»  *  1 


l  +  sir  ,  n  ®  1 ,  •  •  «  t  N 
n 


(2-28) 


As  a  consequence  of  whitening,  (2-25)  can  be  written 


Ur  )  ■  -2.n(l  +  SIR  )  + 
n  n 


SIR 


1  ♦  SIR 


(2-29) 


From  (2-23),  the  log -likelihood  ratio  of  the  whitened  filter  responses, 
r,  is 


N  sir  N 

*<£>•£  1---SIR-  *rJ  -  £  tnu  *  siv  •  (J-30) 

n«l  n»l 

Equation  (2-30)  is  the  usual  form  for  a  quadratic  discriminant. 

A 

Given  this  form,  where  the  expected  value*  of  |r  |  depend  upon  n,  an 
exact,  closed  form  expression  for  system  performance  can  be  found,  and 
this  expression  can  be  written  in  terms  of  the  signal-to- 

interference  ratios  at  the  outputs  of  the  N  orthogonal  filters  obtained 
from  an  SIR  maximization  algorithm. 

By  finding  a  set  of  orthogonal  filters  for  SIR  maximization,  we 
have  obtained  the  major  part  of  an  optimum  detector  configuration.  The 
complete  optimum  detector  is  constructed  by  computing  a  weighted  sum 
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of  magnitude- squared  (square-envelope  detected)  filter  outputs, 

as  in  (2-30).  The  resulting  detector  configuration  is  illustrated  in  Figure  2-1. 

2.3  Detection  Performance 

In  order  to  conform  to  the  notation  in  Van  Trees  [14],  let 

sir  =  A  (2-31) 

n  n 


and 


X 


-1/2 

in 


A 

r 

n 


r 

n 


(2-32) 


as  in  (2-26).  In  this  case 


E  {  |rn|2  |HV)  ■  Xn  ♦  l  (2-33) 

and 

E  {  |  r  j2  |  H  }  »  l  ,  n  *  1,  2#  ...»  N  .  (2-34) 

n  0 

2 

Each  square-envelope  detected  filter  output  |r  |  is  the  sum  of  the 
squares  of  two  uncorrelated  Gaussian  random  variables,  xn  =  Re{rn>  and 
yn  s  lm{rn>  .  The  power  is  assumed  to  be  split  equally  between  these 
two  variables,  so  that  they  have  equal  probability  distributions 


***„»  -  (2ro„2) 1/2  ex‘’['\2/(20..2)] 

p(yn)  -(  2«Jn!)  1/2  a«p[-yn2/(20n;!|] 


(2-35) 

(2-36) 
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Figure  2-1.  Implementation  of  the  Log -Likelihood  Function 
from  Orthogonal  SIR  solutions 


where 


o 

n 


2 


(X  +  i )  / 2  when  H.  is  true 
n  i 


1/2  when  Hq  is  true 


It  follows  that  [21] 


(2-37) 

(2-38) 


piUJ2)  -  p(*n2  +  yn2) 

-  (20  V1  exp(- 1  r  1 2  /  (2a  2)]  .  (2-39) 

n  n  n 

The  corresponding  characteristic  function  is 

!t  1  -  j(Xn  +  l)wj~1  when  H1  is  true  (2-40) 

(1  -  jwj’1  when  Hq  is  true  .  (2-41) 

The  data  dependent  part  of  the  log -likelihood  ratio  is 


Ur) 

m 


£ 


n«l 


♦  1 


(2-42) 


In  order  to  evaluate  detection  performance,  we  want  to  find  the  probability 
distribution  of  2.(r), 
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(2-43) 


PU) 


N 

n 

n  =  l 


C  .F. 


of 


X 

n 

X  +  1 
n 


2 


where  the  double  arrow  indicates  a  Fourier  transform  relation. 


If  the  characteristic  function  of  (rnP  is  <t>n(w),  then  the  C.F. 
of  a(rn|^  is  <J>n(ctw).  It  follows  that 


(l  -  j  X  w]”1  when  H.  is  true 
n  l 


c-F-°f  rrr  I'J2  ■ 

n 


1  '3  X 


X  -I  -l 

-5r-l 

n  J 


(2-44) 


when  Hq  is  true  (2-45) 


and 


N 

n  u  -  j  x  w] 

n«l  n 


-1 


given  H 


1 


(2-46) 


pU> 


N  . 

II  (l  -  )  X'  (Dl  given  H« 
n»l 


(2-47) 


where 


X'  -  X  /(X  +  1) 

n  no 


(2-48) 


The  products  in  (2-46)  and  (2-47)  can  be  represented  as  sums 
by  using  a  Heaviside  expansion,  which  is  extensively  exploited  in  network 
theory  ( 25] .  For  N  distinct  values  of  XR , 
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N 

n 

n  =  l 


(1 


3  An  u) 


E 


a . 

_ 1 _ 

1  -  j  X .  iu 
i 


where 


a . 


x 


1  -  j  X,  a) 


N 

It  (1  - 
n«l 


j  X  a> ) 
n 


w  ■  - j/X^ 


N 

*  n  a 

n-1 

n*i 


W 


Substituting  (2-49)  into  (2-46)  yields 


V  ■  E  \  /  e~ui  d“ 

i  ■  1  -« 

N 

"  X)  exp(-4/Xi) 

i-1 

Similarly, 

N 

p  (  2,  |  Hq  )  -  £  exp(-4/X’) 

i  °  1 


(2-49) 


(2-50) 


(2-51) 


(2-52) 
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where 


N 

a!  =  H  (1 
1  n  =  1 

n/<i 


x;/xi> 


-l 


(2-53) 


Probabilities  of  detection  and  false  alarm  can  now  be  computed  for 
any  threshold  setting,  y.  The  detection  probability  is 


OO 

pd  »  f  pUIV  <n 


N  ® 

°  a  i  J'  exp  (-  (l/X^)  ]  dU/V) 

i»l  Y 

N 

a  expt-Y/^)  (2-54) 

and  the  false  alarm  probability  is 


P f  m  f 

Y 

N 

"  £  al  eKpl-YU^D/XJ  (2-55) 

i  *  1 
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where 


N 

a .  *  H  ( 1 
1  n  =  l 
n*i 


X  /\.) 
n  1 


(2-56) 


and 


N 

n 

n»l 

n^i 


1 


X  /U  +  111’1 
n  n 

xi/txi  +  1)1 


(2-57) 


X.  -  SIR.  ,  (2-58) 

i  i 

the  signal-to-interference  ratio  at  the  output  of  the  iu  filter  f^(x). 

Similar  results  are  found  in  (241 . 


Equations  (2- 54) -(2- 58)  constitute  a  closed  form  expression  for  the 
performance  of  a  detector  for  Gaussian  signals  in  Gaussian  noise,  for  the 
case  in  which  the  signal  covariance  matrix  has  unequal  diagonal  elements. 
If  some  of  the  diagonal  elements  are  equal ,  a  generalized  version  of  the 
technique  in  (2-50)  can  be  used  (251.  The  expression  depends  upon  the 
eigenvalues  of  the  matrix  C(x,y)  in  (2-18).  The  eigenvalues  are  the  same 
as  the  signal-to-interference  ratios  at  the  outputs  of  a  set  of  orthogonal 
filters,  if  the  filters  are  designed  to  maximize  these  ratios. 


2.4  Summary  of  Section  2 

Maximization  of  signal-to-interference  -atio  for  signal -filter  design 
is  generally  used  as  a  way  to  obtain  usable,  albeit  suboptimum,  results. 
A  straightforward  extension  of  the  SIR  algorithm,  however,  yields  not 
only  a  "beat"  signal-filter  pair,  but  a  set  of  additional  orthogonal  filters. 
The  resulting  filter  set  is  apparently  the  same  as  one  would  obtain  with 
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linear  discriminant  analysis  or  with  a  Karhunen-Loeve  transformation  of 
whitened  data.  The  similarity  to  linear  (Kuilback)  discriminants  is  not 
surprising  when  one  considers  that  the  functional  Q(f)  in  (2-5)  is  very 
similar  to  the  Kullback  divergence  [18,19] 

J(f)  =E[«rTf)jH1j  -  E^(rTf)|H0j  (2-59) 

T  T 

where  £(r  f)  is  the  log -likelihood  ratio  when  the  quantity  r  f  is  taken 

as  data. 

The  maximum  SIR  filter  set,  which  also  implements  a  simultaneous 
diagonalization  of  signal  and  interference  covariance  matrices,  can  be  used 
for  optimum  (quadratic  discriminant)  detection.  It  is  only  necessary  to 
form  a  weighted  sum  of  the  squared  envelopes  of  the  filter  outputs.  The 
performance  of  the  resulting  detector  can  be  written  as  an  exact,  closed 
form  expression  that  depends  upon  the  signal-to-interference  ratios  at 
the  filter  outputs. 

3.0  A  POLAR 'METRIC  SCATTERING  FUNCTION 

3.1  Dependence  of  SIR  on  Target  and  Clutter  Scattering  Functions 


From  (2-1),  SIR  depends  upon  the  expected  magnitude -squared 
filter  output 


E 


f*(x)  r(x)  dx 


f*(x)  r^x)  dx 


(3-1) 


where  f^(x)  is  the  filter  that  processes  the  output  of  the  vertically 
polarized  antenna  r^x),  and  f2(x)  is  the  filter  for  the  horizontally 
polarized  antenna  output,  r2(x). 
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The  received  vertically  and  horizontally  polarized  echoes,  r^Cx)  and  r2(x), 
are  related  to  the  transmitted  signal  components  u^x)  and  u2(x)  by  (5] 

2 

r^x)  =  2  /  (x  "  t/2,  t)  uj(x  ‘T)  dT  (3-2) 

j=l  -« 

<*»■« 

where  <t ,  t  )  is  the  impulse  response  of  a  time  varying  random  filter  at 
time  t  when  the  impulse  is  applied  at  time  zero.  In  terms  of  distributed 
radar  reflectors,  (t  *  t/2,  t)  is  the  reflectivity  of  a  scattering  element 
with  delay  t,  measured  at  the  time  of  reflection  t  -  t/2.  The  subscripts 
of  b^  (t,t);  1  =  1,  2;  j  =  1,  2,  imply  that  there  are  actually  four  time 
varying  weights  in  a  tapped  delay  line  model  of  the  filter  at  delay  t,  i.e., 
a  2  x  2  scattering  matrix  that  depends  upon  time  t  and  delay  t.  These 
weights  describe  .he  backscatter  with  polarization  i  for  an  incident  signal 
component  with  polarization  j. 

Substituting  (3-2)  into  (3-1)  yields 


jl  J  £*U)  dx  r 


2  « 


=  e  i  fi*(x)  £  f  vx ' T/2,  x)  uj(x  ' T>  dt 

(  jal  -• 


dx 


-si «,<*-*>  uscy  ■  t,) 

i=l  j=l  m=l  n=l  -» 


b„  (X  -  t/2.  t)  b^n(y  -  T'/2,  t‘)J  dx  dy  dT  dT’  . 


ij 


I 


(3-3) 
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The  Fourier  transform  in  time  of  the  time -varying  tap  weight  at  delay  x 
can  be  used  to  assess  the  Doppler  spread  induced  by  the  time  variation. 
In  order  to  express  the  receiver  output  in  terms  of  Doppler  spread,  one 
can  use  the  definitions 


or 


b^U ,  x) 


bj.(t  -  x/2,  x)  exp(-j2iT<J>t)  dt 


by(x  -  x/2,  x) 


b.j($,  x)  exp(j2ir4>x)  d$ 


(3-4) 


(3-5) 


The  integral  in  (3-4)  is  performed  on  a  sample  function  of  the  time-varying 
tap  weight  b^(t  -  x/2,  x).  The  expected  product  of  two  such  sample 
functions  is 


E 


bjj(x  -  x/2,  x)  b*n(y  -  VI 2,  x') 


00  / 

■  ff  E 

-00  \ 

expCj2xr (4>x  -  $'y)l  d4*d<^’  (3-6) 

This  equation  can  be  simplified  by  assuming  that  tap  weights  at  different 
delays  are  statistically  uncorrelated  and  that  the  temporal  variation  of  the 
weights  is  wide-sense  stationary,  so  that  the  expression  in  (3-6)  is  a 
function  only  of  the  difference  (x  -  x/2)  -  (y  -  x'/2).  This  wide-sense 
stationary,  uncorrelated  scatterer  (WSSUS)  assumption  implies  that 
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I, 


.  .  ) 

=  E  jbyU.T)  b*n(4,T)j  6(4)  -  $')  6( t  -  T») 


and 


E  jb^Cx  -t/2,  t)  bmn(y  -  t'/2,  t») 


/*  (  I 

=  J  E  jbij(4>.T)  bmn^,T)  j  exp[j2ir4>  (x  -  y)] 


(3-7) 


d4>  6( t  -  t') 


(3-8) 


In  (3-8)  the  expectation  is  a  function  of  specific  4>  and  t  values: 


Eiy*,T)  =  E{b„b:M  i  *,t 


mn 


ij  mn 


|  Pb<«t>.T  ) 


(3-9) 


where  pb(4>»x)  is  the  probability  that  the  specified  range  and  Doppler 
values  will  actually  occur. 
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Substituting  (3-9)  into  (3-8)  and  (3-8)  into  (3-3)  yields 


/  »  v  2  2  2  2  ®  oo 

E  |  /  f(x)  r  (x)  dx  |!  // 

v  -«  )  i=l  j=l  m=l  n=l  -»-« 


f  * (x)  u.(x  -  t)  exp(j27T(J)X)  dx 


\l 
U. 

E  |bij  bmn  I***’  Tj  Pb(<J>,x)  d(J)dT 


fm(y)  u*(y  -  t)  exp  ( -J2iT4>y)  dy 


2  2  2  2 


EE  EE// 

i=l  j=l  m=l  n=l  1  J 


Xf  (t,$)  dTd$ 
ran 


where 


Xf  ..  (  t,  4> ) 


f^(t)  u*(t  -  t)  exp  (-j2Tt<t»t)  dt 


is  the  narrowband  cross -ambiguity  function  of  the  reference  function 
and  the  signal  u£(t)  and 

SW*''C)  5  E{btibmn  I  <.t|pb<*.t) 
is  a  polariraetric  version  of  the  target  scattering  function. 


(3-10) 


(3-11) 

:k(t) 

(3-12) 
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3.2  Comparison  with  Less  General  Problem  Formulations 


Fundamental  results  are  embodied  in  (3-10)-(3-12) .  Eq.  (3-10)  shows 
the  dependence  of  SIR  on  the  scattering  function  of  target  or  clutter  and 
upon  signal-filter  design  as  manifested  by  various  cross -ambiguity  functions. 
Eq.  (3-12)  defines  the  scattering  function  in  terms  of  expected  values  of 
products  of  Fourier  transformed  tap  weight  variations,  as  defined  in  (3-4). 

If  i  =  j  =  m  =  n  =  1,  then  S^j(<|>,t)  is  the  power  spectrum  of 
the  time  variation  of  the  tap  weight  at  delay  t,  for  the  vertically 
polarized  channel.  If  only  the  term  i  =  j  =  m=  n  =  lis  considered  in 
(3-10),  then  we  obtain  the  usual  expression  for  SIR  for  the  non-polarimetric 
case  [11].  Polarimetric  processing  introduces  fifteen  additional  terms  into 
the  SIR  expression,  and  significant  improvements  in  SIR  should  occur 
if  there  are  any  polarization -sensitive  differences  between  target  and 
clutter. 

Another  very  general  aspect  of  the  problem  formulation  is  the 
inclusion  of  Doppler  spread  in  cross  polarization  terras,  e.g.,  i  i  j  and/or 
m  i  n  in  (3-10).  These  cases  account  for  twelve  of  the  sixteen  terms  in 
(3-10),  i.e.,  all  the  terms  except  for  i,  ],  m,  n  equal  to  (1111),  (2222), 
(1122),  and  (2211).  Even  if  the  nu«.«uer  of  cross  polarization  terms  is 
effectively  halved  by  the  realistic  assumption  that  by  =  bjj  ,  there  are  still 
six  Doppler  sensitive  cross  polarization  terms  which  may  be  different  for 
target  and  clutter.  If  any  of  these  terms  is  different  for  target  and 
clutter,  the  SIR  maximization  technique  will  exploit  it.  The  most  obvious 

o 

cross  polarization  term  is  E  { |  b12(4> , x )  |  }  ,  the  power  spectrum  of  the 
time  variation  of  b^Ct.x).  Cross -spectral  terms  such  as  E  {  b12U,x) 
b 22(41. x)  )  may  also  be  important. 

It  is  difficult  to  visualize  specific  advantages  of  the  general 
formulation,  to  conceptualize  problems,  and  to  obtain  physical  insight 
into  the  solutions.  In  order  to  facilitate  this  process,  it  is  helpful  to 
use  a  dipole  model  of  the  tap  weights  by(t,t).  This  model  is  discussed 
in  the  next  section. 
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4.0  DIPOLE  MODELS  FOR  ELEMENTARY  SCATTERERS 


4.1  Single  Dipoles 

Instead  of  the  usual  collection  of  point  seatterers,  've  assume  that 
we  have  a  collection  of  dipoles  or  thin  wires.  Let  u^  and  u2  the  vertical 
and  horizontal  components  of  the  transmitted  signal,  and  let  i^  and  ic  be  the 
along-length  and  cross-length  currents  induced  on  the  dipole  by  the  signal. 
For  an  ideal  dipole,  we  assume  that  the  cross-length  current  is  negligible. 
For  a  dipole  that  is  tilted  0  radians  from  vertical,  it  follows  that 


lC 


cos  9 


sin  0 


u- 


(4-1) 


When  the  dipole  re- radiates  energy,  the  vertical  and  horizontal 
components  of  the  echo,  r^  and  ?2>  are  given  by 


cos 

0 

« 

-sin  0 

1l' 

a 

sin 

0 

cos  0 

fc 

j 

m 
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CQS^  0 

■ 

m 

cos  0  sin  0 

U1 

cos  0  sin  0 

sin^  0 

u2 

m 

(4-2) 


where  o  is  the  radar  reflectivity  of  the  dipole.  The  radar  cross  section 
is  f2. 
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4.2  Collections  of  Dipoles 


Point  scatterers  can  also  be  replaced  by  collections  of  dipoles. 
Suppose ,  for  example ,  that  each  uncorrelated  scatterer  is  a  pair  of  dipoles 
that  are  separated  slightly  in  range.  The  scattering  matrix  for  each 
uncorrelated  scatterer  is  then 


S  =  a. 


cos^  0, 


cos  6  ^  sin  8^ 


*  2 


cos  8^  sin  0^  sin  8^ 


+  a-, 


cos^  8. 


cos  02  sin  82 


cos  02  sin  82  sin  02 


.->«<¥•)  .  (4-3) 


where  the- first  dipole  is  oriented  0^  degrees  from  vertical,  the  second 
dipole  is  oriented  02  degrees  from  vertical,  and  the  two  dipoles  are 
separated  by  AR  meters.  In  Eq.  (4-3),  f  is  the  signal  frequency  and 
c  is  speed  of  light. 


In  the  case  of  a  vertical  dipole  in  front  of  a  horizontal  dipole  with 
AR  equal  to  a  quarter  wavelength  and  =  o2,  we  have 


(4-4) 
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The  scattering  matrix  in  (4-4)  induces  a  polarization  reversal  of  a  circularly 
polarized  wave  with  vertical  component  =  A  cos  ( 2ir ft)  and  horizontal 
component  u2  =  A  sin  (27tft). 

If  a  distributed  reflector  is  modeled  as  a  collection  of  single  dipoles 
as  in  (4-2),  then  no  frequency  dependent  phase  shift  is  introduced  into 
the  echo,  since  the  WSSUS  assumption  eliminates  interaction  between 
different  scatterers.  To  obtain  frequency  dependent  phase  shifts,  a 
distributed  reflector  must  be  constructed  from  elementary  reflectors  that 
are  themselves  collections  of  dipoles,  as  in  (4-3).  This  observation  could 
lead  to  a  useful  discriminant  if  the  target  size  is  known,  i.e.,  if  AR  is 

specified  in  (4-3),  and  if  clutter  scatterers  can  be  modelled  as  in  (4-2). 

More  specific  target  models  can  be  obtained  by  considering  specific 
structures  known  to  exist  on  particular  targets,  and  modeling  these 
structures  in  terms  of  measured  scattering  matrices  or  as  combinations  of 
dipoles  and  planar  point  targets.  Planar  point  targets  are  discussed  in 
the  next  subsection. 

4.3  Planar  Point  Target  Models 

Another  type  of  elementary  scatterer  is  the  usual  planar  point 
target  or  perfect  mirror.  Since  the  currents  induced  on  such  a  reflector 
are  parallel  to  the  applied  field,  there  are  no  polarization  shifts.  The 
scattering  matrix  is  oj,,  where  I, is  the  identity  matrix;  b^  =  b22  =  0 
and  b12  =  b2l  =  0.  Thus, 

!E(o2)  if  i  =  }  and  m  =  n 

(4-5) 

0  otherwise 
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As  in  Eq.  (4-3),  one  can  combine  a  dipole  and  a  planar  point 
target  to  model  various  physical  effects.  If  a  planar  point  target  with 
reflectivity  a  is  placed  a  quarter  wavelength  behind  a  vertical  dipole 
with  reflectivity  2a,  then 


1 


S  =  2a 


0 


1 


=  a 


0 


(4-6) 


which  is  the  same  as  the  scattering  matrix  for  circular  polarization  reversal 
in  (4-4). 


A  horizontal  dipole  and  a  ground  plane  with  variable  distance  AR(t) 
between  dipole  and  ground  yield 
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°1 


0 

0 


2AR(t) 

c 


1 


(4-7) 


The  scattering  matrix  model  in  (4-7)  could  be  a  useful  representation  of  a 
low-flying  cruise  missile  against  a  background  of  ground  clutter,  if  both 
missile  and  clutter  are  within  the  same  resolution  cell  [  26] . 


4.  4  Polarimetric  Scattering  Functions  for  Dipole  Scatterers 

Equation  (4-2)  indicates  that  a  dipole  scatterer.  can  be  represented 
by  its  reflectivity  a  and  its  tilt  6  relative  to  vertical.  If  both  of  these 
quantities  are  dependent  upon  Doppler  and  range,  then  the  polarimetric 
scattering  function  in  (3-12)  is 
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Sijmn(*,T)  =  E  jbij  [ 


o(<f, t),  6(4>,t) 


bmn  6(4),  t) 

I  4>»x 

Pb(4>,  T) 

{  Y  r  I 

=  j  j  J  b..(o,C )  b*n(o,e)  p(a,0  |  4>,  t)  dad6 


Pb(4.  t) 


(4-8) 


From  (4-2), 

2 

bj^(a,9)  =  a  cos  9 

bl2(o,0)  =  b  21(o,9)  «  a  cos  9  sin  9 
o 

b22(a,0)  =  a  sin  8 


(4-9) 

(4-10) 

(4-11) 


Assuming  that  a  and  6  are  statistically  independent  random  variables, 
we  have 


p(o,9  |  4>.t)  =p(o  j  <M  >  P(8  I  <M)  (4-12) 

and  Sijmn(4>»T)  is  proportional  to 

00 

E(c2  I  <t>. t)  =  f  o2  p(0  I .t )  do  .  (4-13) 

A  second  simplifying  assumption  is  that  p(8  |  <f>  ,t  )  is  uniformly  and 
symmetrically  distributed  between  ±a(4>,x).  In  this  case, 
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a(<J>,  t) 


Siimn(^*T)  =  E(°2  i  *’T)  Pb(^’T)  f  byO)  b*n(6)  de 


2a(<M) 


a(<j>,i) 


_1_ 

4a 


+  sin  (2a)  *  22n£i°L 


1 

w 


ifi  =  j  =  m=  n=  l 


Y~  -  sin  (2a)  +  8i^(4ot^- 


s  E  (cr  |  4> » x )  pb(4»,x) 


1 

W 


ifi  =  j=  m=  n  =  2 
sin  (4a) 


a 


(4-14) 


,  if  i  4  j  and  m  4  n 


or  if  i  =  j  *  m  =  n 


0,  if  i  »  j  and  m  4  n  or  If 


i  4  j  and  tn  =  n 


where  a  =  a(<fc, t).  The  expressions  on  the  right  hand  side  of  (4-14)  are 
obtained  by  integrating  by  from  (4- 9) -(4-11)  over  9  from  0  =  -a  to  0  =  a, 


Each  tap  weight  in  a  tapped  deiey  line  model  of  an  extended  target 
is  envisioned  as  a  moving  dipole.  In  many  physical  situations,  one  would 
expect  the  Doppler  spread  induced  by  dipole  rotation  to  be  correlated  with 
the  amount  of  rotation,  in  which  case  a($,T)  will  increase  monotonically 
with  $  <  In  any  case ,  the  dipole  model  allows  one  to  conceptualize  the 
difference  between  polarimetric  extended  target  models  and  non-polarimetric 
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models.  We  have  also  obtained  a  specific  expression  for  the  polarimetric 
scattering  function  Such  a  specific  expression  is  very  useful 

for  synthesis  and  analysis  of  SIR  optimization  programs. 

5.0  IMPLEMENTATION  OF  A  COMPUTER  PROGRAM  FOR 
POLARIMETRIC  SIR  OPTIMIZATION 


In  order  to  set  up  a  computer  program  for  SIR  maximization ,  it  is 
useful  to  express  the  signal-to-interference  ratio  in  a  form  that  can  be 
easily  evaluated  by  a  computer  with  user- supplied  target  and  clutter  data. 

To  find  the  best  signal  and  filter  functions  with  a  computer,  it  is  necessary 
to  parameterize  these  functions  in  terms  of  (say)  time  samples  or  frequency 
domain  samples  (Fourier  coefficients).  The  optimization  problem  is  formulated 
by  expressing  SIR  in  terms  of  these  parameters. 


5.1  Signal  and  Filter  Parameters 


The  vertical  component  of  the  transmitted  signal  is  u^tt)  and  the 
horizontal  component  is  Ugft).  Each  component  is  represented  as  a 
weighted  sum  of  complex  orthonormal  basis  functions  (8^(0;  k  -  0,  l,  . ...  K). 
The  vertical  and  horizontal  signal  components  are  thus 


ut(t) 


uikek«)  .  i  - 1.2 


Kk1  V0 


i  a  1,2 


(5-1) 


where  |u^|  is  the  magnitude  of  the  kth  expansion  coefficient  u^  and 
uik  corresPontling  phase  parameter. 
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Similarly,  the  filter  components  are  written 


K  K 

V0  =  Jj,  fik  ek(t)  =  kXj,  |fJ  exp(ivik)  V°  •  (5'2) 

i  =  1,2  . 

The  basis  functions  of  particular  interest  are  either  sinusoids 

QjJt)  =  T”1/2  exp(j2irkt/T),  (5-3) 

or  sin(t)/t  functions  from  sampling  theory 

0^(0  =  B*^  sinc{tTB[t  -  k/B)}  .  (5-4) 

In  (5-3),  T  is  the  signal  or  filter  duration.  In  (5-4),  B  is  the  system 
bandwidth.  Fourier  series  representations  using  the  components  in  (5-3) 
are  desirable  if  the  receiver  already  incorporates  DFT  operations,  as  in 
coherent  pulse  Doppler  or  synthetic  aperture  radars.  Sine  functions  are 
desirable  if  the  radar  uses  a  coded  waveform  and  the  receiver  can 
implement  a  matched  filter  for  such  a  waveform  in  the  time  domain. 

5.2  General  Expressions  for  Expected  Filter  Output  Power  in 

fteaponse  to  target .  Clutter T  of  Noise 

The  expected  filter  output  power  in  response  to  the  target  echo  is 
given  by  (3-10)  In  terms  of  cross  ambiguity  functions  and  a  polarimetric 
scattering  function.  The  ambiguity  functions  are  defined  in  (3-11),  and 
the  scattering  function  is  given  in  (3-12)  and  (4-8).  Similar  expressions 
apply  to  the  filter  output  power  in  response  to  the  clutter  echo,  except 
that  (3-12)  becomes 
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(5-5) 


Sijmn  (+•'')  =  Ej'ij  Cl*’T|pc<4,'T) 

where  and  cmR  are  elements  of  the  clutter  polarization  scattering 
matrix . 


The  expected  filter  output  power  in  response  to  noise  n(t)  is 
E { |  yi*(t)n(t)  dt| 2  =  (Nq/2)  dt 


2 

=  (N0/2)  2  Xf£  (0,0)  .  (5-6) 

i=l  i  i 

5.3  Expressions  in  Terms  of  Expansion  Coefficients  and 
Basis  Functions 


This  section  contains  mathematical  details  that  can  be  ignored 
without  the  loss  of  much  understanding.  If  the  reader  is  not  interested 
in  such  details,  he  can  turn  to  Section  5.4. 


5.3.1  Ambiguity  Functions 


For  orthonormal  basis  functions  as  in  Eqs.  (5-1)  and  (5-2),  we  have 


X*M  (t,<J>) 

A.U. 

l  J 


i  l  “jkfU  / 


k=0  11=0 


For  sinusoidal  (Fourier)  basis  functions  as  in  (5-3),  we  have 
T 


J 9k(t  •'T)9ic(t)  e"j27T4,t  dt  s 


ej2rrkx/T  -j 27r<f>T  . 


<r _  t  - 

2fljU-k  -  $T)le 


1] 


ej((2rkT/T)  -  nri  sincW) 

i  +  t(k-A)/*rr 
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4f 


(5-7) 


(5-8) 


and  for  sinc(t)  basis  functions  we  have  (using*  Parseval's  Theorem) 


_  QQ 

Xfl  (T,*)  =  /  0*(t  -  T)e.(t)  e'j2lt<|,t  dt 

Vk 


/ 


B/2 

-B/2 


(ek(f)  e"i2lTfT]*  +  <»)  df 


r 

-B/2 


®i(f  +  ^  ®k(f)  eJ2irfT  df  ‘ 


(5-9) 


where  6^.(f )  is  the  Fourier  transform  of  9k(t),  i.e 

•1/2 


0k(f)  = 


B  exp(-j2fffk/B)  ,  -B/2  <  f  <  B/2 

0  ,  otherwise 


(5-10) 


It  follows  that,  for  sine  basis  functions, 


XA  a  ( ^.0 )  -  e  0/B)  gjncfiKk-i,  +  B*d) 

Vk 


(5-11) 


Expressions  for  the  ambiguity  functions  in  (3-11)  are  thus  as 
follows.  For  sinusoids  on  a  time  interval  { 0 » T ] , 

u,  f  ei2itkt/T 

X£  u(‘t*(^)  s  e  ^  ^  sinc(ir<$T)  E  E*  I'-V'nk-fl f^]  .  (5-12) 
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and  for  sine  basis  functions  with  frequency  domain  support  on  [-B/2,  B/2], 


Xf  -  zz  ujf.j  e  j27T^/B  sinc(ir(k-4  +  Bt)]  .  (5-13) 

iuj  k  SL  a 


5.3.2  Expected  Filter  Output  Power  in  Response  to  Time -Varying, 
Distributed  Target  or  Clutter 


Substitution  of (t,<|>)  from  (5-12)  or  (5-13)  into  (3-10)  yields 
iu! 

expressions  for  the  expected  filter  output  power.  These  expressions  contain 
the  integral 


(5-14) 


Further  analysis  of  this  integral  is  possible  if  it  is  assumed  that  p  (t  ,  ) 
can  be  represented  by  a  two  dimensional  histogram,  i.e.,  p(t  ,<J> )  is 
constant  over  histogram  intervals  of  size  At,  A4>.  The  integral  in  (5-14) 
is  then  a  sum  (over  the  indices  r  and  s)  of  simpler  integrals  of  the  form 


6+A4/2  t  +At/2 
s  /•  r  - 


Pb(V*.,E!b  /  /  X*f  <’■*>  *(„<''*)  dT  d« 

*A*n  t/m/2  ‘  i  m  “ 


(5-15) 


The  double  integral  in  (5-15)  can  be  evaluated  by  using  (5-12) 
for  sinusoidal  basis  functions  or  (5-13)  for  sine  functions.  For  sinusoids, 


(5-16) 


<t>  +A<i>  12  t  +At  12 
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(4>T)^sinc^(7T<J>T) 

[(j)T  +lk-l")]  UT  T  (p-q)l 


d<$> 


and 


<j>g- A  <|>  /  2 


(5-17) 


Tr+Ar/2 


l2  = 


f  e-i2*<k-p>T/T  dT  _ 


T.-AT/2 


The  integrals  in  (5-17)  and  (5-18)  are  evaluated  in  Appendix  A. 


(5-18) 


For  sine  functions,  we  have  the  same  form  as  in  (5-18)  except  that 


y*-At/2 


I,  = 


J r  sinc(  TT(k-£  +  Bt)]  sinc[Tr(p-q  +  Bt)jdT 


(5-19) 


y-At/2 
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m-j  «uau-i 


and 


<j>  +A<J>/2 

1 2=  J  e-j2*U-q)$/B  .  (g 

4>s-A4>/2 

The  integrals  in  (5-19)  and  (5-20)  are  also  evaluated  in  Appendix  A. 

From  (3-10),  the  expected  filter  output  power  in  response  to  the 
target  echo  rT(t)  is 


E{|/f*(t)rT(t)dt  |2} 


X>b<v*s>  £  E(bijb;jw 


£  u„.f  *  u*.f. 


r,s 


l.j.m,  n 


k, H.p.q 


a  u  a 

jk  iS,  np  mq 


I1(k,«.,p,q.Tr,(|)s)I2(k,2,,p,q,Tr,<j>s) 


(5-21) 


A  similar  expression  yields  the  clutter  response,  provided  that 

EIcijcmn  lYV  is  used  instead  of  E[bijbmnlTr,<,,s1,  and  pb(VV  18 

replaced  by  PC<V  V# 

From  (5-6)  and  (5-7),  the  expected  filter  output  power  in  response 
to  noise  is 


E(i/y  (On(t)dt  | 2 }  = 


<Nq/2> 


2  K 

!C  £ 

i=l  k=0  lK 


(5-22) 
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5.4  Signal-to- Interference  Ratio  Optimization 


Using  Eqs.  (5-21)  and  (5-22),  the  signal- to- interference  ratio  can 
be  written 


SIR  = 


E{  1  target  response  } 


•y 

E{|clutter  response  |  /  +  E{|noise  response  I*'} 


(5-23) 


where 


E  ( |  target  response  |2}  Pb(VV  X*  E(bijbmn  IVV 


r,s 


S  V*l  u*„p£mq  ’  <5-2« 


k,£,p,q 


E{|dutter  response]2}  =  £pc<V*s>  £  E(cij4n 


r,s 


i.j.m.n 


£  Vil  “np  mq  tl(k-t'P><»'V'>s)I2(k'1'>’’‘1' W 


(5-25) 


k, i.p.q 


2  K 


E  { | noise  response  p)  =  (NQ/2)  £  £  lfik 

i=i  k-0 


(5-26) 


The  optimization  problem  is  to  find  the  signal  and  filter  parameters 
( luilc  I  *  uik’  l^ik  I  *  vik*  1  -  i»  2;  k  a  0 . K)  such  that  SIR  is  maximized. 

From  the  above  equations,  it  is  straightforward  to  obtain  the 
derivatives  of  SIR  with  respect  to  each  unknown  coefficient,  and  a  gradient- 
search  type  of  optimization  program  can  be  used  to  maximize  SIR.  In 
Section  2,  however,  it  was  shown  that  a  likelihood  ratio  test  is  not 
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generally  implemented  with  a  single  linear  filter;  a  sequence  of  orthogonal 
filters  must  be  used.  Each  filter  in  this  sequence  should  maximize  SIR 
for  the  signal  specified  by  the  first  signal-filter  pair,  and  each  filter 
function  should  be  orthogonal  to  the  previously  determined  filter  functions. 

Such  a  sequence  of  orthogonal  filters  can  generally  be  obtained  by 
adding  constraints  to  the  SIR  expression  via  Lagrange  multipliers  [19,27], 
An  easier  technique,  however,  is  to  exploit  a  property  of  the  simplex 
method  for  maximizing  a  function  [  16] .  The  starting  or  trial  solutions 
in  the  simplex  algorithm  define  a  subspace  containing  the  "optimum” 
solution  found  by  the  algorithm.  To  constrain  the  search  to  functions 
that  are  orthogonal  to  a  set  of  previous  solutions,  th->  initial  points 
(simplex  vertices)  can  deliberately  be  made  orthogonal  to  the  previous 
solutions,  and  no  further  modifications  of  the  algorithm  are  necessary. 

A  disadvantage  of  the  simplex  technique  as  it  presently  exists  in  the 
literature  is  that  it  functions  without  gradient  information.  Such  informa¬ 
tion  is  available  for  SIR  and  would  presumably  speed  up  the  search  for  a 
maximum  if  it  were  used.  A  summary  of  the  discussion  in  [16]  about  the 
simplex  method  is  given  below. 

A  simplex  is  the  convex  hull  of  n+1  points  in  Rn:  A  triangle  in 
2  3 

R  ,  a  tetrahedron  in  R  ,  etc.  If  we  want  to  maxmiize  a  function  of  n 
variables,  we  evaluate  the  function  at  each  of  n+l  vertices  which  specify 

t. 

an  initial  simplex.  For  one  of  these  points,  x  ,  the  function  f(x)  is 

0 

largest,  and  for  another  point,  x  ,  it  is  smallest.  The  object  at  each 
step  is  to  replace  x  ,  the  vertex  of  the  current  simplex  with  the  lowest 
function  value,  by  a  new  and  better  point. 

A  tentative  direction  for  the  new  point  is  obtained  by  drawing  a 
line  from  x^  through  the  mean  of  all  the  other  points,  computed  by 
excluding  x2'.  The  resulting  point  is  computed  from  the  "reflection" 
operation 
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(5-27) 


xr  =  x  +  a(x  -  x2)  , 

where  a  is  a  positive  constant  (a  value  of  1  is  recommended)  and 


*  =  5  S.*i  .  -***  •  (5-28) 

1=0 

If  f(xh)  <  f(xr),  i.e.,  the  reflection  step  has  generated  a  new  maximum, 
then  we  take  an  "expansion”  step  in  the  same  direction  by  computing 

xe  =  x  +  y(xr  -  x)  ,  (5-29) 

where  y  >  1  is  a  given  constant  (a  value  of  2  is  recommended).  If 
f(xe)  >  f(xr),  then  xe  replaces  xr.  If,  however,  f(xe)  <  f(xr),  then 

a  - 

the  expansion  step  failed  and  x  is  replaced  by  x  to  form  a  new  simplex. 

If  the  reflection  step  results  in  a  new  point  that  no  longer  has  the 
smallest  function  value  but  is  also  not  the  largest  value,  then  x  is 
replaced  by  xr  without  implementing  an  expansion  step. 

Another  possibility  is  that  the  original  reflection  step  fails  in  the 
sense  that  the  point  xr  still  has  the  lowest  function  value,  i.e., 

f(xr)  <  min  (f(x.)  ,  Xj  ^  .  (5-30) 

i 

In  this  case,  xr  would  just  be  a  new  version  of  the  least  desirable  point. 

If  (5-30)  is  true,  then  a  "contraction"  step  is  used: 

xc  =  x  +  Six*1  -  x)  ,  (5-31) 
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where  0  <  $  <  1  is  the  contraction  coefficient  (a  value  of  8  =  1/2  is 
recommended)  and  x^  is  identical  with  either  xr  or  x^,  depending  upon 
which  point  yields  the  largest  function  value,  i.e., 


o ' 


f(xc  )  =  max  {f(x  ),  £(x  )} 


(5-32) 


Vi  c 

If  x  =  x  ,  then  x  establishes  a  new  vertex  in  the  direction  opposite 
to  xr,  as  we  see  by  comparing  (5-27)  and  (5-31).  If  =  xr,  then 
we  proceed  in  the  direction  of  xr  -  x  from  the  point  x  ,  but  our  step 
size  is  0  <  8  <  1  rather  than  y  >  1  as  in  (5-29). 

If  f(xc)  <  f(x 4  ),  then  we  still  need  to  find  a  new  vertex  that  is 
not  the  least  desirable  point,  i.e.,  we  still  have  not  reversed  the  inequality 
(5-30).  In  this  case,  all  the  original  simplex  vertices  are  moved  toward 
the  point  with  the  largest  function  value,  i.e,  the  whole  simplex  is 
concentrated  near  the  best  point.  The  new  simplex  vertices  are: 


A  1  ,  h  , 
*i  =  xi  *  2  (x  ’  xi) 


i  =  0,  . . . ,  n 


(5-33) 


A  suggested  termination  criterion  is  based  on  the  observation  that 
as  a  result  of  (5-33),  all  the  vertices  of  the  simplex  near  a  maximum  will 
move  close  together  and  close  to  x.  Thus,  if 


n 


rh  £  - f(x)l2 

i=0 


(5-34) 


is  sufficiently  small,  then  the  algorithm  should  terminate. 


An  important  property  of  the  simplex  method  is  its  dependence 
upon  the  initial  vertices  chosen  for  the  first  simplex.  For  example,  if 


these  vertices  do  not  completely  span  Rn,  then  the  solution  will  be  found 
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in  a  subspace  of  Rn.  Although  this  property  may  be  viewed  as  an 
annoyance,  it  is,  in  fact,  very  useful  for  finding  solutions  that  are 
orthogonal  to  a  given  function,  e.g.,  a  previously  determined  filter 
function.  The  importance  of  this  observation  has  been  pointed  out  with 
respect  to  filters  for  maximizing  SIR  and  the  optimum  (likelihood  ratio 
test)  receiver  structure. 

5.5  Iterative  Optimization  of  Signal  and  Filter  Vectors 

In  order  to  find  the  best  signal-filter  pair  for  SIR  optimization, 
it  is  nearly  always  necessary  to  iteratively  improve  one  set  of  parameters 
with  the  others  held  fixed.  We  have  four  functions  to  optimize:  The 
vertically  polarized  signal  u^(t),  the  horizontally  polarized  signal  u2(t), 
the  vertically  polarized  filter  function  f^(t),  and  the  horizontally  polarized 
filter  function  f2(t).  Each  of  these  functions  has  been  represented  in 
terms  of  two  parameter  sets:  the  magnitudes  of  the  expansion  coefficients 

{|uikl,  I  u2kl*  I  f ik  I »  |f2kl  *  k  =  °*  1 . K)  and  the  corresponding 

phase  parameters  (ulk»  U2k>  v2^;  k  =  0,  1,  ....  K}.  There  are 

thus  eight  vectors,  each  with  K  +  1  terms,  that  specify  the  signal  and 
filter : 


lujJ,  ]u2 1 ,  IfjJi  !£_2 1  *  Hi*  ^2*  Hi*  H.2 

The  vertical  and  horizontal  signal  component  magnitudes,  |  uj  and  |  u2 1 » 
are  linked  by  an  energy  constraint: 

2  K 

Signal  Energy  =  £  |  u,.  |  2  =  1  (5-35) 

i=l  k=0 

The  coefficients  |  u^|  and  |  u2  I  are  scaled  in  order  to  satisfy  (5-35)  before 
they  are  used  to  evaluate  SIR  at  each  stage  of  the  simplex  algorithm.  The 
search  for  an  optimum  signal  is  thus  constrained  to  the  space  of  unit 
energy  functions. 
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The  link  between  1 |  and  |u2|  in  (5-35)  implies  that  a 
straightforward  optimization  technique  should  find  all  the  signal  magnitude 
coefficients  (vertical  and  horizontal)  with  other  coefficients  (filter 
magnitude,  signal  phase,  and  filter  phase)  held  fixed.  The  same  is  true 
of  the  filter  magnitude  coefficients  if  one  uses  the  constraint 

2  K  „ 

Filter  Energy  =  E  E  !*«,!  =  1  •  (5-36) 

i=l  k=0  ** 

It  can  be  argued  that  a  constraint  on  filter  energy  is  unnecessary  because 
SIR  is  insensitive  to  multiplication  of  f(t)  by  a  nonzero  scalar  constant. 

For  ease  in  interpreting  the  results  as  well  as  for  faster  convergence, 
however,  (5-36)  has  been  implemented  as  well  as  (5-35).  The  resulting 
simplex  algorithm  automatically  energy  normalizes  all  signal  or  filter  magnitude 
coefficients  before  evaluating  SIR,  and  the  method  iteratively  optimizes 
signal  magnitude,  filter  magnitude,  signal  phase,  and  filter  phase  with 
all  other  components  held  fixed. 

6.0  TESTING  THE  SIR  MAXIMIZATION  ALGORITHM:  A  SIMPLE  EXAMPLE 

This  section  describes  a  simple  test  of  a  computer  optimization 
program  for  maximization  of  signal-to- interference  ratio  (SIR).  The  test 
demonstrates  the  basic  input /output  parameters  of  the  program.  The 
results  are  reasonable  from  an  analytical  viewpoint,  indicating  that  the 
algorithm  is  working  properly. 

The  test  is  incomplete  in  two  important  respects  which  will  be 
addressed  further  on  in  the  report.  First,  there  is  no  polarization- 
dependent  random  Doppler  or  frequency  spread  difference  between  target 
and  clutter  (or,  alternatively,  no  Doppler- dependent  polarization  effects). 
Second,  only  the  first  SIR  filter  is  found,  whereas  we  have  shown  that 
additional  filters  can  be  used  to  implement  a  Bayes  optimum  detector 
(likelihood  ratio  test)  for  discriminating  signal  from  interference. 
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6.1  Description  of  the  Simple  Test  Problem 


The  polarimetric  scattering  function  for  dipoles  that  are  uniformly 
distributed  in  vertical  tilt  between  ±a(<f> ,  t)  is  given  by  (4-14). 

As  a  simple  test  case,  let  the  target  dipoles  be  nearly  vertical 
(-1°  <  a  <  1°)  and  let  the  clutter  dipoles  be  oriented  randomly 
(-90°  <  a  <  90°).  In  this  case,  if  E ( a2 1 » x )  =  1,  the  target  scattering 
function  is 


(<t> .  t  ) 


Pb(d>*T) 


if  i  =  j  =  m  =  n  =  1 


otherwise  , 


(6-1) 


and  the  polarimetric  scattering  function  of  the  clutter  is 


if  i  =  j  =  m  =  n  s  1  or  2 

if  i  *  j  and  m  *  n 

if  i  =  j  *  m  =  n 
otherwise  ,  (6-2) 

where  pb(<j>,t)  is  the  distribution  of  the  target  in  Doppler  and  range,  and 

p  ( <j> , x )  is  the  corresponding  clutter  distribution.  In  order  to  accentuate 
c 

discrimination  on  the  basis  of  polarization,  it  will  be  assumed  that  pb(4>,T) 

equals  p  (<J>,x).  This  assumption  eliminates  features  that  are  generally  used 
c 

for  non-polarimetric  SIR  maximization. 


(3/8)  pc(<t>,x) 


(1/8)  P0(<t,T) 


or 
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6.2  Analytical  Investigation  of  the  Simple  Test  Case 
The  signal -to-interference  ratio  is 

_  E[  [filter  response  l2|target  echo  alone] _ 

E[  [fltr.  resp.  | 2 (clutter]  +  E[  jfltr.  resp.  | 2 (noise] 

Substituting  ( 6- 1)  into  ( 3- 10) 

0° 

E{  |fltr.  resp.  |2  (target)  =  ff  Pb($,t)|xf  u  C  )  (at  d<{> 

11 

Substituting  (6-2)  into  (3-10) 


E { {fltr.  resp.  | 2 (clutter }  =  ~  a  PC(<M)  E(t,4>)  dt  d-J» 


where 


£<  t,$) 


“  I  !xf1u1<T’l),)  t2  *  I  l*f2u2<I'»,lZ 

*  5  litflu2<,’*,lZ  *  I  l*f2u1(,-*)|i 

♦  g  (2  RelXf  u  (T.*>  X,  u 

••  1  b  ft 


(6-3) 


(6-4) 


(6-5) 


47 


+  |  {2  Re[x!  (t.$)  Xf  u  (*.$)]} 
I1U2  X2U1 


=  illxf  .ul2+  lxf  f23 

1  1  2  2 


*5tlW  V2'*  +  !V2  +  V/1  ' 


(6-6) 


Finally,  the  noise  response  is 

9  N 

E[  |fltr.  resp.  |  | noise )  =  -y  ( Xf  <■  (0.0)  +  Xf  f  (0.0)] 

“  xri  X2X2 

S  N  12  .  (6-7) 

o 

In  the  computer  algorithm,  NQ  is  set  equal  to  0.02  in  order  to  allow  the 
SIR  to  be  dominated  by  clutter.  The  filter  energy  in  (6-7)  is  constrained 
to  be  unity,  as  is  the  signal  energy  (the  sum  of  the  squares  of  the 
vertical  and  horizontal  signal  components). 

Even  with  all  our  simplifications,  it  is  still  difficult  to  analytically 
obtain  a  signal* filter  pair  to  maximize  SIR.  It  is  possible,  however,  to 
suggest  two  different  solutions  on  the  basis  of  physical  and  mathematical 
insight. 


From  a  physical  viewpoint,  we  observe  that  the  target  dipoles  are 
nearly  vertical,  while  the  clutter  dipoles  are  uniformly  distributed  in 
orientation.  From  this  viewpoint,  we  would  expect  all  signal  and  filter 
energy  to  be  concentrated  in  the  vertical  direction,  while  the  horizontal 
components  should  have  zero  energy.  Our  physical-insight  solution  is 
then 
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(6-8) 


u2(t)  =  f2(t)  =  0 


which  implies  that 


//  ph(<{>,T)ixf  u  d<J> 

A  -  00  *'  1  ‘  1 

SIR  = - - - LA - _ - 

3/8  // lxf  u  rdx d<i> +  *01 
-*>  11 


if 


2.67 


(6-9) 


Pb(4>.T)  =  Pc(4>.t)  .  (6-10) 

From  a  mathematical  viewpoint,  a  possible  solution  is  obtained  by 
noting  that  the  clutter  response  in  (6-5)  is  minimized  if  fj(t),  f2(t), 
u^(t)  and  u2(t)  are  chosen  so  as  to  minimize  the  right  hand  side  of  (6-6). 
To  make  the  expected  clutter  response  small,  let 


fiui  =  - 


and 


Xf2u2 


X.'. 


nu2 


xf  u 
r2ul 


(6-11) 

(6-12) 


such  that  the  two  last  terms  in  (6-6)  are  both  minimized.  In  order  for 
(6-11)  and  (6-12)  to  hold  simultaneously,  we  can  let 

u1(t)  =  -  u2(t)  ,  (6-13) 

f^t)  =  f2(t)  .  (6-14) 
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A  simple  sign  reversal  or  180°  phase  shift  is  thus  a  tentative  solution.  In 
this  case, 


SIR 


//  Pb(<t>»T)  lx;  ~  (t,<*> )  [2  dx  d<{> 

_ -j* _ 11 _ 

CO 

kf[  Pe(*.t>l IxJj^l2  +  IXf2~2l2]  d*  dr*  .01 


(6-15) 


where  (6-13)  and  (6-14)  imply  that 


(t.*j I2  =  lx;  (t.*> I2 

1U1  J2U2 


Substituting  (6-16)  into  (6-15),  we  have 


SIR  s  2 


(6-16) 


(6-17) 


where  pb(<J>,x)  =  pc(4>.x)  as  in  (6-10). 


Physical  insight  and  the  form  of  the  clutter  response  in  (6-5) 
and  (6-6)  have  suggested  two  possible  solutions  to  the  SIR  maximization 
problem.  The  first  solution  yields  the  larger  SIR,  but  this  solution  is 
only  expected  to  be  optimum  when  the  target  dipoles  are  vertical  and  all  the 
clutter  dipoles  are  horizontal.  In  the  given  problem,  clutter  dipole 
orientations  are  not  all  horizontal,  but  are  uniformly  distributed  with 
respect  to  tilt  angle. 


6.3  Computer  Solution  for  the  Test  Case 


A  computer  solution  for  the  test  case  has  been  generated  with  the 
algorithm  described  in  Section  5.  From  Eq.  (3-12)  the  target  scattering 
function  is 


slW*,T)  =  E{bij  bmnl*'T)  V*lT) 


(6-18) 


where 


/  b..(o,e)  b*  (o,9) 
o  =  0  mn 


pb( a, 0 (<J> , t)  do  d6  .  'b_1 

Expression  (6-19)  is  the  description  of  the  target  that  is  employed  in  the 
SIR  maximization  algorithm.  In  addition  to  (6-19),  one  must  also  specify 
the  distribution  Pb(4i,T)  as  in  (6-18).  Similar  expressions  are  used  for 
the  clutter. 

The  input  to  the  computer  program  was  thus  E{by  as 

given  by  (4-14)  for  the  target,  with  a  =  it/ 180°  radians.  For  the  clutter, 
E{cjj  c*n|$,t}  was  also  given  by  (4-14),  but  with  a  =  tt/2.  The 
distributions  Pb(4»,t)  and  pc(4»,t)  were  both  uniform  on  the  t.4>  plane, 
extending  from  -T  to  T  in  range  and  from  -4-rr/T  to  4ir/T  in  Doppler, 
where  T  is  the  duration  of  the  signal. 

The  iterations  that  were  performed  by  the  computer  are  described 
and  documented  in  Appendix  B.  Two  runs  were  made.  In  the  first  run, 
the  magnitudes  of  signal  and  filter  Fourier  coefficients  were  adjusted 
first,  with  the  phases  held  fixed.  The  phases  were  then  adjusted  with 
magnitudes  held  fixed,  etc.  In  the  second  run,  the  phases  were  adjusted 


51 


first.  The  results  are  in  terms  of  the  Fourier  coefficients  of  the  signal 

and  filter  functions,  i.e.,  the  parameters  lfik'  *  yik’  Vik;  1  =  1,2; 

k  =  0,  ....  4}  in  the  expansions  (5-1)  and  (5-2).  For  the  first  run, 
the  vertically  polarized  signal  and  filter  components  of  the  solution  are 


u10  =  0.317  e]7T 

f1Q  =  0.307  e^° 

u..  =  0.411  e^ 
li 

fn  =  0.516 

u12  =  0.599  e^ 

f12  =  0.555  e]° 

ul3  =  0.359  e1’77 

f13  =  0.353  e]0 

u..  =  0.281  ei7T 

14 

f14  =  0.141  e^° 

horizontally  polarized  components  are 

u20  =  0.H3  e*0 

f20  ■  °-u9  e)0 

u21  =  0.139  e^° 

f21  =  0.333  e*° 

u22  =  0.293  e3° 

10 

f22  =  0.191  e J 

u23  =  0.138  e5° 

JO 

f23  =  0.128  6 

u24  =  0.141  ej0 

f24  =  0.097  e5° 

(6-20) 


(6-21) 
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The  resulting  SIR  is  2.697,  which  is  better  than  either  SIR  in  (6-9)  or 
sTr  in  (6-17).  The  solution,  however,  seems  to  be  a  compromise  between 
Eq.  (6-8)  and  Eqs.  (6-13)  and  (6-14).  The  horizontally  polarized  signal 
and  filter  functions  in  (6-21)  have  less  energy  than  the  vertically  polarized 
functions  in  (6-20),  but  the  horizontal  energy  has  not  been  driven  to 
zero  as  in  (6-8).  The  vertical  signal  has  been  multiplied  by  exp(jiT)  =  -1 
as  in  (6-13),  but  (6-11)  is  obviously  violated. 

The  second  run,  which  starts  by  optimzing  coefficient  phases  rather 
than  magnitudes,  yields  somewhat  different  coefficient  values  (see  Appendix  B). 
Nevertheless,  the  same  multiplication  of  u^(t)  by  -1  and  a  similar  imbalance 
of  energy  in  favor  of  the  vertical  components  is  again  observed,  and  the 
SIR  is  2.680. 

7.0  DETECTION  OF  AN  EXTENDED  TARGET  IN  SEA  CLUTTER 

In  this  section,  the  SIR  maximization  algorithm  is  applied  to  a  more 
challenging  problem  with  some  practical  significance.  The  goal  is  to  apply 
the  method  to  a  problem  that  involves  Doppler  dependent  polarization 
modulation .  Such  a  problem  arises  naturally  in  the  context  of  radar  sea 
echo.  The  dipole  modelling  concept  in  Section  4  is  especially  useful  in 
this  context. 

Sea  clutter  can  be  said  to  possess  Doppler  dependent  radar 
reflectivity  and  polarization  properties,  Doppler  spread,  however,  is  not 
really  the  independent  variable;  radar  reflectivity  and  polarization  both 
depend  upon  wind  speed,  and  so  does  Doppler  spread. 

According  to  M.  Skolnik's  Radar  Handbook,  p.  26-14:  "In  calm 
seas  with  little  wind,  the  echo  obtained  with  horizontal  polarization  is 
considerably  less  than  that  with  vertical  polarization.  The  echo  with 
horizontal  polarization  increases  with  increasing  wind  speed  faster  than 
the  increase  with  vertical  polarization,  so  that  with  rough-sea  conditions 
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there  is  less  difference  in  the  magnitude  of  the  echo  from  horizontal  or 
vertical  polarization."  [6] 

According  to  p.  26-30  of  [61,  some  measurements  indicate  that  the 
spectral  width  of  sea  echo  is  approximately  proportional  to  wind  speed. 
Quantities  that  vary  with  wind  speed  might  then  be  said  to  vary  with 
Doppler- induced  frequency  spread.  If  the  sea  is  modelled  by  a  collection 
of  dipoles,  then  the  orientations  and  cross  sections  of  these  dipoles  are 
correlated  with  wind  speed,  but  we  can  say  that  "polarization  and  cross 
section  is  Doppler  dependent"  if  Doppler  spread  (<j>),  rather  than  wind 
speed,  is  viewed  as  the  independent  variable. 

A  calm  sea  can  be  represented  with  dipoles  that  are  randomly 
distributed  over  a  small  interval  in  vertical  angle  (8); 

-<x<*)<6  <a(*)  (7-1) 

where  8  is  the  dipole  tilt  measured  from  vertical  and  a(<t>)  is  the  maximum 
excursion  from  vertical.  As  wind  speed  increases,  8  increases  to  as  much 
as 

-ir/2  <  8  <  tt/2  ,  (7-2) 

yielding  random  polarization  for  high  wind  speed.  The  dipole  cross  section 
(or  the  density  of  dipole  reflectors)  also  increases  with  wind  speed  and  thus 
with  Doppler  spread  (#). 

A  boat  or  ship  can  be  modelled  as  a  collection  of  randomly  oriented 
dipoles  or  as  planar  reflectors  (specular  point  targets).  These  dipoles  or 
specular  glints  can  be  assumed  to  be  uniformly  distributed  over  a  range 
interval  that  is  small  relative  to  the  clutter  extent. 

Translational  motion  of  the  target  causes  its  mean  Doppler  frequency 
to  be  displaced  from  the  clutter  mean,  while  pitching  and  rolling  cause  the 
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target  echo  to  exhibit  some  Doppler  spread  about  the  mean.  The  radar 
cross  sections  of  the  dipoles  or  specular  glints  comprising  the  target  are 
independent  of  wind  speed  and  thus  of  <t>. 


It  will  be  of  interest  to  determine  the  effect  of  the  dipole  model  vs. 
the  planar  reflector  target  model  upon  receiver  design  and  performance. 

If  a  significant  difference  exists,  then  it  should  be  possible  to  discriminate 
a  target  composed  of  planar  reflectors  from  a  chaff  cloud  composed  of 
randomly  oriented  dipoles.  For  planar  point  scatterers,  the  scattering 
matrix  is  where  £  is  the  identity  matrix;  b^  =  b22  =  0  while  b12  *  b21 
It  follows  that,  for  planar  point  targets, 


The  models  used  here  portray  the  qualitative  description  of  sea 
clutter  in  [6]  with  the  simplest  possible  functions:  Uniform  probability 
distributions  and  linear  dependence  on  Doppler  spread. 


which  corresponds  to  a  uniform  distribution  of  dipole  reflectivity  a 
between  zero  and  10  dB. 


The  distribution  of  the  target  on  the  range-Doppler  plane  completes 
the  target  scattering  function  description  in  (4-14).  Let 


V*-T)  *[£  *"»  (SiTT1)]  [4  rect  (j^)]  . 


(7-6) 


where 


(  1  .  -1  1  x  <  1 

rect(x)  = 

*  0  ,  otherwise 

(7-7) 


The  target  acatterer  distribution  has  a  Doppler  spread  of  2ir/T  on  either 
aide  of  a  mean  Doppler  frequency.  The  mean  Doppler  shift  is  equal  to 
2tt/T,  corresponding  to  translational  motion.  The  signal  duration  is  T, 
so  that  2tr/T  denotes  a  Doppler  resolution  cell  or  bin  width.  The  delay 
spread  of  the  target  is  2T/5  seconds,  or  two  range  resolution  cells.  A 
range  resolution  cell  for  a  single  frequency  component  T  seconds  long  is 
%T  seconds.  Five  frequency  components  yield  five  times  the  bandwidth,  and 
the  range  resolution  cell  is  thus  approximately  T/5  seconds.  A  top  view  of 
p^(<Ji.T),  looking  down  on  the  range-Doppler  plane,  is  shown  in  Figure  7-1. 
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4  (DOPPLER  SPREAD) 


7.1.2  Specular  Glint  Target 

A  collection  of  specular  (planar  target)  glints  is  specified  by 

2 

E{b^  b*n  in  (7-3)  and  by  pb(x»<i>).  It  can  be  assumed  that  E(ab) 

equals  10/3  as  in  (7-5)  and  that  pb(t,<i>)  is  the  same  as  in  (7-6)  and 
Figure  1. 


7.1.3  Sea  Clutter  at  Low  Wind  Speed  (Calm  Sea) 

Let  the  maximum  tilt  of  the  clutter  dipoles  be  given  by  a  linear 
function  of  Doppler  magnitude: 


10ir  80tt  l<t> 

“  180  TSC  (Wf) 


(7-8) 


The  average  clutter  dipole  cross  section  is  defined  as 


E{oc2|«)  »$[ 


1  + 


(7-9) 


which  is  another  linear  function  of  Doppler  magnitude.  The  sea's  polarization 
spread  and  reflectivity  are  thus  assumed  to  be  linearly  dependent  upon 
Doppler  spread,  which  is  itself  monotone  increasing  with  wind  speed.  If 
the  true  variation  of  maximum  dipole  tilt  ac  with  wind  speed  W  is  ac  =  g(W), 
and  if  Doppler  spread  is  related  to  wind  speed  by  |$i  =  h(W),  then  our 
simple  model  assumes  that  g(h  l(|$|)l  is  linear. 
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Plots  of  otc( 4> )  and  E{o  (<!>}  are  shown  in  Figures  7-2  and  7-3.  For  a 
calm  sea,  these  plots  are  only  relevant  for  |<j>|  <  2tt/T  ,  since  it  is  assumed  that 


Pc(*’T)  =  47rect 


(7-10) 


i.e.,  a  calm  sea  has  no  Doppler  spread  beyond  |<j>|  =  2tt/T.  The  area  of 
the  T,<j>  plane  covered  by  p  (<j>,t)  is  shown  in  Figure  7-4. 

The  assumption  that  p_< <f> , x )  is  uniformly  distributed  on  (-T ,T) 
in  the  delay  (t)  direction  is  equivalent  to  an  assumption  of  uniformly  dis¬ 
tributed  clutter  for  all  ranges.  This  equivalence  follows  from  the  fact  that 
the  ambiguity  function  is  nonzero  only  over  a  delay  interval  between  -T  and 
T  when  the  target  is  hypothesized  to  be  at  delay  zero. 


e(oc2!$} 


4 -  ROUGH  SEA  - » 


Figure  7-2.  Assumed  variation  of  clutter  dipole  cross  section  with 
maximum  clutter  Doppler  spread 
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Figure  7-3.  Assumed  maximum  effective  dipole  tilt  relative  to  vertical 
(for  a  distributed  dipole  model  of  sea  clutter)  as  a 
function  of  maximum  clutter  Doppler  spread. 


4>  (DOPPLER  SPREAD) 


Figure  7-4.  Assumed  delay -Doppler  distribution  of  sea  clutter  for 
low  wind  speed  (calm  sea) 
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7.1.4  Sea  Clutter  at  High  Wind  Speed  (Rough  Sea) 


For  a  distributed  dipole  model  of  sea  clutter,  it  has  been  assumed 
that  the  maximum  tilt  of  any  dipole  is  ac(<f>) ,  and  that  the  tilt  is  uniformly 
distributed  between  -a  (<t»)  and  a (4)*  To  be  consistent  with  the 

C  v 

observations  in  [6],  the  maximum  tilt  should  increase  monotonically  with  4 
until,  at  maximum  wind  speed  and  Doppler  spread,  a  (4  )  equals  tt/2 . 

v  UlaX 

At  maximum  wind  speed,  vertically  and  horizontally  polarized  dipoles  are 
equally  likely,  and  there  is  no  preferred  polarization  for  minimizing 
sea  echo. 

A  simple  linear  dependence  of  ac(4)  on  |<t>|  can  accomplish  the 

above  objectives.  This  depeiidence  is  the  same  as  in  (7-8),  provided 
that 

IW  ■  8*/T  •  (?-u> 

The  resulting  maximum  tilt  as  a  function  of  4  is  shown  in  Figure  7-3. 

The  experimentally  observed  increase  in  sea  clutter  cross  section 
with  increasing  wind  speed  or  Doppler  is  modelled  as  in  (7-9).  For 
Doppler  shifts  limited  as  in  (7-11),  the  expected  clutter  cross  section 
varies  between  0.33  and  3.33,  as  shown  in  Figure  7-2. 

The  assumed  distribution  of  sea  clutter  in  delay  and  Doppler  for 
a  rough  sea  is: 

Pc(4.t)  =  rect  [ff  rect  (t)]  ’  <7*l2) 

so  that  the  maximum  Doppler  spread  is  given  by  (7-11).  The  area  covered 
by  this  distribution  is  shown  in  Figure  7-5. 
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Figure  7-5,  Assumed  delay -Doppler  distribution  of  sea  clutter  for 
high  wind  speed  (rough  sea) 
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The  SIR  maximization  program  has  been  applied  to  signal- filter 
design  for  the  situations  given  in  Table  1.  In  most  cases,  only  one 
filter  has  been  associated  with  the  optimum  signal.  It  was  shown  in  Section  2, 
however,  that  multiple  orthogonal  filters  should  be  obtained  in  order  to 
implement  a  Bayes  optimum  detector  as  well  as  one  with  maximum  signal  to 
interference  ratio.  The  additional  filters  are  easy  to  find  by  always  using 
the  signal  that  yields  largest  SIR  with  the  first  filter,  i.e.,  the  signal 
associated  with  the  "best"  signal-filter  pair.  Additional  filter  functions, 
orthogonal  to  the  ones  found  previously,  are  then  determined  such  that 
SIR  is  maximized  with  the  given  signal.  Since  a  simplex  algorithm  is  used, 
orthogonalization  is  accomplished  by  using  initial  "guesses"  or  simplex 
vertices  that  are  all  orthogonal  to  previous  solutions.  The  solution  from 
such  an  initial  simplex  is  in  the  subspace  defined  by  the  "guesses,"  and  is 
thus  orthogonal  to  previously  derived  filter  functions  for  SIR  maximization. 
Multiple  orthogonal  filters  were  found  for  two  cases:  the  specular  glint 
(planar)  distributed  target  in  rough  sea  clutter  and  the  same  target  type 
in  uniformly  distributed  dipoles  (chaff).  The  latter  problem  is  described 
in  Section  7.5.  The  last  line  of  Table  1  refers  to  a  polarimetric  inter¬ 
ference  canceller  to  be  discussed  in  Section  7.6.  This  device  computes 
the  weighted  difference  between  the  outputs  of  a  signal  detector  with 
maximum  SIR  and  a  clutter  detector  with  minimum  SIR. 


TABLE  1 


COMPUTER  EXPERIMENTS 


Target  Model 

Clutter  Model 

Receiver  Configuration 

Random  dipole 

Calm  sea 

One  filter 

(H  and  V  components) 

Random  dipole 

Rough  sea 

One  filter 

(H  and  V  components) 

Specular  glint 

Calm  sea 

One  filter 

(H  and  V  components) 

Specular  glint 

Rough  sea 

One  filter 

(H  and  V  components) 

Specular  glint 

Rough  sea 

Many  orthogonal  filters 

Specular  glint 

Random  dipole 
(chaff) 

One  Filter 

(H  and  V  components) 

Specular  glint 

Random  dipole 
(chaff) 

Many  orthogonal  filters 

Specular  glint 

Random  dipole 
(chaff) 

Polarimetric  interference 
canceller 

t 


7.3 


Results  for  the  R.sndom  Dipole  Target  Model  in  Calm  and 
Rough  Seas 


If  the  target  is  a  collection  of  randomly  oriented  dipoles  and  the 
clutter  is  composed  of  nearly  vertical  dipoles,  the  problem  is  very  similar 
to  the  simple  one  that  was  analyzed  in  Section  6.  The  solution  for  sea 
clutter  would  then  be  a  system  with  energy  concentrated  in  the  horizontal 
signal  and  filter  channels,  with  one  signal  component  phase  shifted  by  180°. 


In  the  sea  clutter  problem,  a  difference  in  range-Doppler  clutter 
distributions  has  been  introduced.  The  target  has  a  measurable  average 
Doppler  displacement  relative  to  the  clutter,  and  the  range  extent  of  the 
target  is  restricted.  These  two  changes  should  ideally  result  in  signals 
with  large  time-bandwidth  product.  If  the  best  strategy  is  to  exploit  only 
Doppler  resolution,  however,  then  the  signal  energy  will  become  concentrated 
at  a  single  Doppler  component,  and  the  filter  energy  will  do  the  same. 

In  summary,  a  reasonable  solution  would  involve  horizontally 
polarized  signal  and  filter  functions,  with  signal  or  filter  phase  shifted  by 
180°  and  with  energy  concentrated  at  one  frequency  to  exploit  Doppler 
displacement  between  target  and  clutter. 


A  similar  solution  has  been  obtained  by  the  SIR  maximization 
algorithm.  For  a  calm  sea,  the  vertically  polarized  signal  Fourier  components 
{Uiq.  uu....  ,uu)  and  vertically  polarized  filter  components  {f^,  •  ♦  • » f  14 ^ 

are  relatively  small: 


10 

'  0. 

003 

expi  jr) 

f10 

11 

=  0. 

037 

exp'  ]r) 

fll 

12 

=  0. 

007 

exp(jr-) 

f12 

13 

=  0. 

015 

e.\p(  }•») 

f13 

14 

_  °- 

005 

exp(jr) 

f14 

0.009  exp(jO) 

0.002  exp(jO) 

0.003  exp(jO)  (7-14) 

0.005  exp(jO) 

0.003  exp(jO) 
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The  horizontally  polarized  signal  and  filter  components  for  a  calm  sea  are 


u20  *  °*222  exp(jO)  f 2q  =  0.454  exp(jO) 

U21  =  ^*602  exp(jO)  f^  =  0.755  exp(jO) 

u22  =  0.630  exp(jO)  f22  =  0.469  exp(jO) 

u23  =  0.407  exp(jO)  f23  =  0.056  exp(jO) 

u2,  =  0.152  exp(jO)  f24  =  0.003  exp(jO) 

and  the  resulting  signal-to-interference  ratio  for  a  calm  sea  is 

SIR  =  51.54  . 

In  Eq.  (7-15),  there  appears  to  be  a  mismatch  between  signal  and 
filter,  such  that  the  filter  pass  band  is  approximately  a  frequency  shifted 
version  of  the  transmitted  signal  spectrum.  The  shift  is  approximately  one 
frequency  component  downward.  From  Fig.  7-1,  the  average  target  Doppler 
shift  is  exactly  one  frequency  resolution  cell,  2tt/T.  A  positive  <J>- value 
in  Fig.  7-1  corresponds  to  a  target  with  range  increasing  with  time,  which 
means  that  the  echo  is  shifted  downward  in  frequency  by  an  average 
Doppler  displacement  corresponding  to  one  frequency  component.  A  filter 
matched  to  the  expected  echo  will  then  have  a  transfer  function  that  is  a 
downward -shifted  version  of  the  signal  spectrum,  as  observed  in  (7-15). 

The  simplex  SIR  maximization  algorithm  used  three  iterations  for  adjustment 
of  each  set  of  coefficients  (signal  phase,  filter  phase,  signal  magnitude, 
filter  magnitude).  The  total  run  time  for  computer  design  of  the 
polarimetric  radar  was  about  six  hours,  but  the  program  used  disk 
storage  for  many  computed  variables.  Subsequent  examples  were  computed 
using  large  RAM  arrays,  and  run  time  was  reduced  to  approximately 
two  hours. 

For  a  rough  sea,  slightly  more  energy  is  relegated  to  the  vertically 
polarized  signal  and  filter  components,  as  one  would  expect  from  Figure  7-3. 
For  the  calm  sea  case,  part  of  the  target  Doppler  spread  is  outside  the 


(7-15) 


(7-16) 
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Doppler  band  of  the  clutter,  as  shown  by  Figures  7-1  and  7-4.  For  the 
rough  sea,  the  target  Doppler  spread  is  wholly  immersed  in  the  clutter, 
as  shown  by  Figures  7-1  and  7-5.  This  difference,  along  with  larger 
average  clutter  reflectivity,  leads  to  a  smaller  maximum  SIR  for  rough  seas. 

The  vertically  polarized  signal  and  filter  components  for  the 
rough  sea  model,  as  obtained  by  the  SIR  optimization  algorithm,  are: 


u^0  =  0.019  exp(jTr) 
Uji  =  0.122  exp(jTr) 
u12  =  0.061  exp(jir) 
u13  =  ^  .017  exp(jrr) 
u14  =  0.021  exp(jrr) 

The  horizontally  polarized  signal  and 

u2q  =  0.203  exp(jO) 
uf.  =  0.758  exp(jO) 
u22  s  0.598  exp(jO) 

u23  =  °-076  exp(jO) 
u24  =  0.030  exp(jO) 


f^Q  =  0.011  exp(jO) 

fj^  =  0.010  exp(jO) 

f^2  =  0.044  exp(jO)  (7-17) 

f13  =  0.002  exp(jO) 

f14  =  0.005  exp(jO) 

filter  components  for  a  rough  sea  are: 

f2Q  =  0.422  exp(jO) 

f21  =  0.809  exp(jO) 

f22  -  0.406  exp(jO)  (7-18) 

f23  »  0.002  exp(jO) 

f24  =  0.003  exp(jO) 


The  above  signal  and  filter  functions  yield  a  signal -to-interference 
ratio  for  the  rough  sea  clutter  model  of 


SIR  =  26.83 


(7-19) 


As  one  would  expect,  this  SIR  is  considerably  less  than  for  a  calm  sea. 

The  signal  and  filter  functions  for  SIR  maximization  for  a  random 
dipole  target  in  sea  clutter  combine  Doppler  resolution  with  accentuation 
of  horizontally  polarized  energy.  As  in  the  simple  test  case  in  Section  6, 
the  vertical  signal  components  are  all  multiplied  by  minus  one  or  exp(jtr). 
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This  result  is  not  generally  exploited  by  sea  surface  search  radars,  and 
the  design  may  yield  improved  performance  in  rough  seas. 

7.4  Results  for  the  Range- Distributed  Planar  Target  (Specu^r  Glint) 
Model  in  Calm  and  Rough  Seas 

The  specular  glint  target  model  is  similar  to  the  classical 
non-polarimetric  concept  of  a  distributed  target  composed  of  randomly 
spaced  "highlights'*  or  points  of  high  reflectivity  that  behave  as  perfect 
mirrors  (planar  reflectors).  These  highlights  are  assumed  to  have  the 

same  cross  section  as  the  corresponding  dipoles  in  the  distributed  dipole 

o 

target  model,  i.e.,  E[o  ]  is  the  same  in  (7-3)  and  (4-14). 

2 

The  same  values  of  E[cr  )  in  (7-3)  and  (4-14),  however,  seem  to 
make  the  planar  target  more  detectable.  For  example,  if  the  target  dipoles 
are  randomly  tilted  such  that  a  equals  ir/2  in  (4-14),  we  have 


Sijmn(<M) 


randomly 

tilted 

dipoles 


*  E(o4)  p($,t) 


as  in  (6-2),  while  from  (7-3), 


3/8  if  i  a  ]  ~  m  =  n  =  1  or  2 

1/8  if  i  *  j  and  m  *  n 

or  if  *  =  j  ^  m  =  n 

0  otherwise 

(7-20) 


Hjran 


(<t>,  t) 


-  E(a  )  p($,T> 


planar 

points 


1  ifi  =  j  =  m=  n  =  lor2 
1  if  i  *  j  t  m  a  n 
0  otherwise 

(7-21) 


If  most  signal  and  filter  energy  is  horizontally  polarized,  then  the 
dominant  term  in  the  polariraetric  scattering  function  is  $2222^,T^’ 
is  8/3  times  larger  for  the  planar  target  model  than  for  uniformly 
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distributed  dipoles. 


Despite  the  difference  between  planar  targets  and  randomly  tilted 
dipoles  in  (7-20),  a  non -random  dipole  orientation  is  consistent  with  (7-21). 
For  a  vertically  oriented  dipole  with  a<<l,  (4-14)  yields 


Sijmn 


(<t>,  T ) 


=  E(a2)  p(«M) 

vert 

dipole 


1  if  i  =  j  =  m  =  n  =  1 
0  otherwise 


(7-22) 


which  is  commensurate  with  the  planar  point  target  in  (7-21).  A  randomly 
oriented  dipole  has  smaller  effective  cross  section  than  the  same  dipole 
with  known  orientation.  The  planar  point  target  cross  section  has  been 
chosen  to  be  equal  to  that  of  a  dipole  with  known  orientation,  while  SIR 
depends  upon  the  effective  cross  section  of  a  randomly  oriented  dipole. 


The  above  observations  imply  that  the  effective  radar  cross  section 
of  a  target  can  be  reduced  by  more  than  3  dB  if  different  reflecting 
surfaces  have  different  polarisations  or  effective  dipole  tilts.  If  all 
reflecting  surfaces  have  the  same,  known  polarization,  then  the  equivalent 
specular  glint  model  must  use  a  higher  glint  cross  section  in  order  to 
represent  the  target. 


For  a  calm  sea,  the  SIR  maximization  algorithm  with  a  distributed 
planar  target  model  yields  vertical  signal  and  filter  coefficients 


IQ3 

0.010 

exp 

(jO) 

ll3 

0.019 

exp 

(jO) 

52 

12 

0.002 

exp 

(jo) 

13" 

0.001 

exp 

(jo) 

14= 

0.014 

exp 

(jo) 

f10=  0.007  exp  (jO) 
f^=  0.001  exp  (jO) 
f12-  0.005  exp  (jo) 
f13=  0.004  exp  (jO) 
fi4=  0.003  exp  (jO) 
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and  horizontal  coefficients 


U20= 

0.229 

exp 

(JO) 

f20= 

0.400 

exp 

(3'0) 

U21= 

0,590 

exp 

(JO) 

f21= 

0.779 

exp 

(j0) 

U22= 

0.606 

exp 

(j0) 

f22= 

0.481 

exp 

(j0) 

U23= 

0.426 

exp 

(]0) 

f23= 

0.041 

exp 

(jo) 

U24= 

0.223 

exp 

(30) 

f24= 

0.001 

exp 

(jo) 

The  best  SIR  for  the  planar  glint  distributed  target  model  in  calm  seas 
is  found  to  be 

SIR  =  137.87  (7-25) 


When  the  planar  glint  target  model  is  used  with  our  model  for  a 
rough  sea,  »ve  obtain  vertical  coefficients 

ulQ  =  0.028  exp(jir) 

Ujj  =  0.004  exp(jfr) 
ul2  a  0*002  exp(jir) 
u^  =  0.004  exp(jfl) 
u14  =  0.009  exp(Jir) 


and  horizontal  coefficients 

u20  =  0.212  exp(jO) 
u2l  3  0.780  exp(}0) 

u22  =  0,578  exP^°* 

u23  =  0.105  exp(jO) 
u24  -  0.017  exp(jO) 

The  corresponding  SIR  for  a  target 
reflectors  in  a  rough  sea  is 

SIR  =  75.43  . 


f1Q  -  0.006  exp(jO) 

«  0.008  exp(jO) 

fl2  =  C.004  exp(jO)  (7-26) 

f12  »  0.003  exp(JO) 
fl4  »  0.025  exp(jO) 


f2Q  a  0.340  exp(jO) 

t'21  3  0.856  exp(jO) 

f22  -  0.388  exp(jO)  (7-27) 

f23  =  0.001  exp(jO) 

f24  a  0.000  exp(jO) 

with  distributed  planar  (point  target) 

(7-28) 
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From  our  interpretation  of  (7-20)  -  (7-22),  we  expect  that  the  SIR 
for  the  planar  glint  target  model  will  be  8/3  or  2.67  larger  than  the 
corresponding  SIR  for  the  distributed,  randomly  oriented  dipole  target 
model.  Indeed,  for  calm  seas,  the  signal-to-interference  ratio  in  (7-25) 
is  2.68  times  larger  than  the  one  in  (7-16),  and  SIR  in  (7-28)  is  2.81  times 
larger  than  SIR  in  (7-19). 


7.5  Results  for  Distributed  Planar  Reflectors  in  a  Background  of 
Randomly  Oriented  Dipoles 


An  interesting  problem  is  to  discriminate  between  the  two  target 

types:  planar  reflectors  and  randomly  oriented  dipoles.  In  this  case, 

2 

E(o  )  for  the  planar  targets  is  set  equal  to  3/8  in  order  to  compensate 
for  the  disparity  between  (7-20)  and  (7-21).  The  range-Doppler 
distribution  of  both  targets  is  the  same,  and  is  shown  in  Figure  7-1.  The 
only  difference  between  the  two  models  is  that 


(<J>,t) 


randomly 

tilted 

dipoles 


=  p<  t) 


1/8  if  i  i  j  and  m  *  n 
or  if  i  =  J  i  m  =  n 


(7-29) 


while 


a  P<$,  t) 

planar 

points 


Sijmn(<M) 


3/8  if  i  =  j  i  m  »  n 
0  if  i  *  j  and  m  *  n 


(7-30) 


The  SIR  maximization  algorithm  must  thus  rely  upon  differences 

»"  ®U22'  S2211'  S1212-  S1221'  S2112’  and  S2121‘  The  ftrSt  tw0  ,erma 
are  larger  for  the  planar  point  target  model,  and  the  last  four  terms 

are  larger  for  the  randomly  tilted  dipoles. 
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The  SIR  maximization  algorithm  for  the  target  discrimination 
problem  yields  vertical  signal  and  filter  coefficients 


u1Q  =  0.316  exp(jO) 
=  0.316  exp(jO) 
u12  =  0.316  exp(jO) 
u13  =  0.316  exp(jO) 
u14  =  0.316  exp(jO) 

and  horizontal  coefficients 

u2Q  =  0.316  exp(JO) 
u21  =  0.316  exp(jO) 
u22  =  0.316  exp(jO) 
u23  =  0.316  exp(jO) 
u24  =  0.316  exp(jO) 


f1Q  =  0.316  exp(jO) 

f^  =  0.316  exp(jO) 

f12  =  0.316  exp(jO)  (7-31) 

f13  =  0.316  exp(jO) 

f14  =  0.316  exp(jO) 


f2(J  =  0.316  exp(jO) 

f21  *  0.316  exp(jO) 

f22  =  0.316  exp(JO)  (7-32) 

f23  =  0.316  exp(jO) 

f24  =  0.316  exp(jO) 


The  corresponding  signal-to-interference  ratio  is 

SIR  =  0,973  .  (7-33) 


The  above  results  can  be  explained  by  substituting  the  polarimetric 
scattering  functions  sijmn(o*T)  for  target  or  clutter  into  (3-10).  For 
distributed  planar  reflectors,  substitution  of  (7-21)  into  (3-10)  yields 


E(  (target  response!  2}  - 


'f2u2 


p(t*$)  dr  d$ 


(7-34) 
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and  for  distributed  randomly  oriented  dipoles 


E{ | clutter  response  |2}  =  |  Jf^  lxf  u  1 2  +  !xf  u  |23 

IX  2  2 

+  ,Xflul  +  V21  *  'Xflu2  +Xf2ul'  ’ 


p(x ,4>  )  dT  d$ 


(7-35) 


as  in  (6-5). 

From  (7-31)  and  (7-32),  the  SIR  maximization  algorithm  has  found 
a  solution  of  the  form 

ux(t)  =  u2(t)  =  fx(t)  =  f2(t)  i  u(t)  .  (7-36) 


Substitution  of  (7-36)  into  both  (7-34)  and  (7-35)  yields 


SIR  3 


3/2  //lxuu(T,$)|2  p(T,<fr)  dr  d 
3/2  //!xuu(t^)|2  p(T,<t>)  dx  dx  +  Nq/2 


(7-37) 


**»  1. 


as  in  (7-33). 

It  can  b@  shown  from  (7-34)  and  (7-35)  that  the  same  SIR  value 
is  obtained  if  all  signal  and  filter  energy  is  concentrated  in  one 
polarization  channel  (either  horizontal  or  vertical).  This  result  seems  to 
imply  that  polarimetric  radar  is  not  superior  to  non -polarimetric  radar  with 
respect  to  detection  of  planar  target  distributions  in  chaff  composed  of 
randomly  oriented  dipoles,  despite  the  differences  between  (7-29)  and 
(7-30).  It  will  be  demonstrated,  however,  that  polarimetric  radar  is 
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superior  if  a  clutter  cancellation  method  is  used.  The  use  of  many  ortho¬ 
gonal  filters,  as  in  Section  2,  has  also  been  investigated. 

The  interference  cancellation  concept  is  discussed  in  vhe  next 
section,  and  the  use  of  multiple  orthogonal  filters  is  discussed  in 
Section  7.7. 

7.6  A  Polarimetric  Interference  Canceller 

The  poor  detection  performance  implied  by  (7-33)  and  (7-37)  can 
be  greatly  improved  by  augmenting  the  receiver  with  a  second  one,  in 
parallel  with  the  first,  which  estimates  the  interference  response  of  the 
first  receiver.  The  estimated  interference  is  then  subtracted  out  or 
cancelled. 

In  order  to  estimate  interference  in  the  possible  presence  of  a 
target  echo,  the  interference  estimator  should  have  relatively  small  response 
to  the  target.  Such  a  processor  will  have  small  SIB,  Alternatively,  it 
should  have  large  SIR  for  the  inverse  problem  where  target  and  clutter 
are  interchanged  and  the  target  echo  is  viewed  as  interference. 

For  planar  reflector  vs.  random  dipole  discrimination,  the  inverse 
problem,  i.e.,  detecting  an  array  of  randomly  tilted  dipoles  while  minimizing 
receiver  response  to  specular  glints,  yields  much  larger  SIR  than  the 
original  problem ,  This  large  SIR  difference  for  the  original  and  inverse 
problems  is  the  key  to  successful  clutter  cancellation . 

The  interference  cancellation  concept  is  illustrated  in  Figure  7-6/ 
Filter  #2  for  interference  estimation  should  Ideally  be  able  to  observe  the 
interference  in  the  absence  of  signal.  In  other  words,  the  interference-to- 
signal  ratio  should  be  maximized  by  filter  #2,  or  the  output  signal-to- 
interference  ratio. 
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sir2  =  s2/i2 


(7-38) 


r<t) 


should  be  minimized.  Alternatively,  the  SIR  for  the  inverse  problem, 
with  clutter  and  target  interchanged,  should  be  maximized  by  filter  #2. 


A  A 

V!2 


Figure  7-6.  An  Interference  canceller  using  filters  for  SIR 

maximization  and  minimization. 
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Figure  7-7  shows  the  polarimeiric  version  of  Figure  7-6.  In  the  figures, 

A  A 

and  I j  are  estimated  or  predicted  interference  levels  at  the  outputs  of 
filter  #1  and  filter  #2,  respectively. 


Figure  7-7.  Pclarimetric  version  of  the  interference  canceller 
in  Figure  7-6. 


In  order  to  accentuate  the  response  of  the  interference  canceller 
to  the  desired  signal,  the  system  should  maximize 


h  -  <WS2  =  S1 


SJU  ' 

spiT. 


(7-39) 


which  implies  that  filter  #1  should  be  designed  to  maximize  SIR,  while 
filter  #2  should  be  designed  to  minimize  it,  as  already  stated. 


The  expected  magnitude -squared  target  response  in  (7-34)  is  zero 
tor  the  second  filter  in  Figure  7-7  if 


u 


(xA) 


1 


*  -X 


4S,'« 


(xA) 

2 


(7-40) 


where 


Uj(t)  *  u2<t)  *  u(t) 


(7-41) 


from  (7-36).  Eq.  (7-40)  is  satisfied  If 
f<2)(t)  =  *f[2)(t)  «  u(t )  . 


In  this  case,  the  SIR  for  the  inverse  problem  (target  and  clutter 
Interchanged)  is 


(SIR) 


inverse 


1/2  //jxuu(^)i2  dt  d$ 


(7-42) 


(7-43> 


For  the  receiver  defined  by  (7-42), 


S2  F  0 


(7-44) 
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ar.d  from  (7-43)  and  (7-37) 


JjSj'l  •  <7-45> 

The  output  of  filter  #2  should  be  multiplied  by  approximately  three  in  order 
to  cancel  the  interference  in  Figures  7-6  and  7-7. 

The  resulting  inteference  canceller  theoretically  allows  the  polari- 
metric  radar  to  "look"  through  chaff.  From  (7-39)  and  (7-44),  the  response 
to  the  target  echo  or  signal  is  unaffected,  while  the  response  to  interference 
is  greatly  reduced. 

Although  the  example  given  above  does  not  exploit  Doppler  dependent 
cross -polarization  effects,  the  general  concept  can  be  applied  to  a  variety 
of  problems. 

7 . 7  Computer  Results  for  Multiple  Orthogonal  Filters 

The  analysis  in  Section  2  has  demonstrated  that  an  SIR  maximi¬ 
zation  result  can  be  extended  to  yield  a  Bayes  optimum  detector  or  likeli¬ 
hood  ratio  test  implementation  by  adding  orthogonal  filters,  as  in  Figure 
2-t.  It  would  appear  that  the  n—  orthogonal  filter  should  increase  SIR 
a8  well  as  yielding  a  better  approximation  to  a  likelihood  ratio  test.  It 
will  be  shown,  however,  that  an  SIR  increase  does  not  occur,  even  though 
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1L 

the  n —  filter  causes  to  increase  more  than  Pp  for  sufficiently  high 
detection  threshold,  Y.  Although  SIR  is  a  useful  measure  of  performance 
for  system  design,  it  can  be  misleading  as  a  detectability  measure. 


7.7.1  SIR  Computation  when  Multiple  Filters  are  Used 

Equations  (2-54)  and  (2-55)  yield  the  improvement  in  performance 

(P„  vs  Pp)  that  is  expected  to  occur  from  the  use  of  multiple  orthogonal 
Dr 

filters.  Another  measure  of  performance  is  the  signal-to-interference  ratio 
itself,  as  measured  at  the  output  of  the  system  in  Figure  2-1.  To  compute 
this  SIR  value,  one  can  exploit  the  observation  that  the  output  of  each 
orthogonal  filter  is  statistically  independent  of  the  output  of  any  other  filter. 
If  the  filter  outputs  are  denoted  r^  r2>  ....  rn,  and  if  the  n—  filter 
output  is  transformed  with  a  function  fn(rn) ,  then  statistical  independence 
implies  that 

e{  s  w  }  *£  E I  w !  •  (7'46) 

n=l  n=l 


For  the  receiver  in  Figure  2-1, 


W  -^a. 


SIR 


n 


In 


1  +  SIR 


n 


(7-47) 
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where 


SIR 

n 


(7-48) 


The  SIR  of  the  whole  filter  bank  in  Figure  2-1  is  the  output  when  only 
signal  is  present,  divided  by  the  output  when  only  interference  is  present. 
When  only  signal  is  present, 


E  {irJ2  I  }  5 


Sn 


(7-49) 


and  for  interference, 


1 1  r^  |  ^  }  interference  |  =  X 


In 


(7-50) 


It  follows  that 


and 


N 

E  |  output  I  signal  |  =  E  j  J  fn(rn>  I  i 


N 

=  E 

n=l 


SIR 


n 


(7-51) 


1  +  SIR 


n 


E  j1  output  J  interference  | 


N 

E 

n=l 


SIR 


n 


(7-52) 


1  +  SIR. 


n 
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Finally,  the  signal-to-interference  ratio  at  the  output  of  the  system  in 
Figure  2-1  is 


SIR  = 

r  N 

E 

SIBn2  1 

r  N 

E 

S,Rn  | 

-  n=l 

l  +  SIRn  J 

.  n=l 

l  +  SIRn  J 

(7-53) 


The  computer  optimization  program  yields  the  optimum  signal,  a 

set  of  optimum  orthogonal  filters,  and  the  SIR  value  at  the  output  of 

each  filter,  1SIR.I  .  The  SIR  values  can  be  used  to  compute  ROC 
t  n  f  n-l 

curves  from  (2-54)  -  (2-55),  and  to  compute  the  overall  SIR  from  (7-53). 


Although  additional  filters  as  in  Figure  2-1  are  needed  to  implement 
a  likelihood  ratio  test,  (7-53)  implies  that  additional  filters  do  not  increase 
SIR.  It  has  been  shown  that 


SIR1  >  SIR2  >  ....  >  SIRn 


(7-54) 


since  the  SIRn  values  correspond  to  solutions  of  a  maximization^ problem 
with  n-1  orthogonality  constraints,  and  because  SIRn  is  the  n-~  largest 
eigenvalue  ofC(x.y)  in  (2-18).  If  SIRn  <  SIRj  for  all  n  >  1.  then  (7-53) 

becomes 


SIR. 


SIR  =  1  +  SIR1 


1  + 


SIR1 
1  +  SIRj 


+ 


(1  +  SIR1  \ 

1  *  SIRn  / 

(1  +  SIR1  \ 


<  SIR1 


(7-55) 
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since  each  term  in  the  numerator  sum  is  equal  to  each  term  in  the 
denominator  sum,  multiplied  by  (SIRn/SIR^>,  which  is  less  than  unity. 

Eq.  (7-55)  states  that  SIR  for  more  than  one  filter  is  always  less  than  SIR 
for  a  single  "best"  filter.  The  use  of  multiple  orthogonal  filters  as  in 
Figure  2-1  cannot  be  justified  in  terms  of  SIR  improvement.  In  fact,  SIR 
decreases,  indicating  a  decrease  in  the  ratio  of  the  mean  receiver  output 
given  signal  to  the  mean  receiver  output  given  interference. 

SIR,  however,  may  not  be  a  good  measure  of  performance.  The 
usual  detectability  index,  for  example,  is 


d 


2 


[E  U(r)  |  Hx)  -  E  U(r)  |  HQ}  ]2 
Var  U< r)  |  HQ} 


SIR 


j  N 

ln?i  1  +  SIR_  (SIRn 


r  N 

I  2 

|_n=l 


r+ 


SIRn 
SIR 


n 


(7-58) 


For  N=l,  i.e.,  for  a  single  filter  pair  f^(t),  d2  equals  KSIR1  -  l)2, 
which  is  monotone  increasing  with  SIR^  if  SIR1  >  1.  The  d  measure 
decreases,  however,  when  additional  filters  have  SIRn  <  1.  Theoretically, 
such  additional  filters  may  still  improve  detectability.  The  problem  is  that 
ROC  performance  depends  strongly  on  the  tails  of  the  distributions  p  ( fi.  j  H 
and  pU|H0)  in  (2-54)  and  (2-55),  while  SIR  and  d2  are  not  very  sensitive 
to  tail  behavior. 


The  above  results  imply  that  although  SIR  is  a  useful  measure  for 
system  design ,  it  should  be  used  with  caution  for  system  evaluation.  For 
system  evaluation,  the  expressions  for  and  Pp  in  (2-54)  and  (2-55) 
should  be  used.  When  the  n—  filter  is  added  in  Fig.  2-1,  PD  increases  by 
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APD  =  an  exP  (-y/SIRn) 


(7-57) 


and 


increases  by 


APp  =  an’  exp  [-Y(SIRn  +  D/SIRJ 
=  an*  exp  (-Y/SIRn>  exp  (-Y) 


(7-58) 


where  Y  is  the  threshold  level  of  the  detector.  It  follows  that  AP^  >APp  if 


a  >  a  '  exp  (-Y) 
n  n 


(7-59) 


or  if 


exp  (T)  >  a  7a„ 
n  n 


n-1 

n 

1  .  SIRi/(SIRi  +  1) 

i=l 

«V  SIRn  +  i) 

n-1 

1 

n 

1  -  (SIR, /SIR  ) 
l  n 

i=l 

n- 1 

=  SIR  El  (SIR,  +  1)  /SIR. 
n  i=l  1 


(7-60) 


The  addition  of  the  n—  filter  will  improve  receiver  performance  if  the 
threshold  v  is  such  that  (7-60)  is  satisfied.  From  Fig.  2-1,  £(r)  is 
always  -  0,  and  a  negative  threshold  value  makes  no  sense.  Therefore, 
it  always  pays  to  add  an  extra  filter  if  the  filter  output  has  a  signal-to- 
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interference  ratio  such  that  the  right  hand  side  of  (7-60)  is  greater  than 
unity .  This  result  follows  from  expressions  for  and  Pp,  and  is  not  in 
agreement  with  SIR  results. 

7.7.2  Further  Results  for  a  Target  Consisting  of  Distributed  Planar 
Reflectors  in  Randomly  Oriented  Dipole  Clutter 

After  finding  the  signal  and  filter  functions  u,(t),  u0(t)  and  f/^(t), 

/  \  \  14  1 

f0v i;(t)  in  Figure  2-1,  the  initial  simplex  is  modified  to  include  only  filters 

£  ( i) 

that  are  orthogonal  to  f  (t).  The  resulting  filter  function  for  maximum 
( 21 

SIR  is  f  (t),  where  the  same  signal  components  u^t),  u2(t)  are  assumed. 

After  finding  u.(t),  u2(t),  f/^tt),  f2(1)(t),  f1<2)(t> ,  and  f2(2)(t), 

another  filter  f^(t)  is  found,  such  that  f^(t)  is  orthogonal  to  both 

f(1)(t)  and  f^2)(t),  The  n—  filter  f(n^(t),  has  signal -to -noise  ratio 

SIR  ,  and  SIR  should  decrease  monotonically  with  n  because  of  the  extra 
n  n 

orthogonality  constraints  that  are  included  as  n  increases.  Computer  results 
for  a  target  consisting  of  a  random  distribution  of  planar  reflectors  and  for 
clutter  consisting  of  randomly  oriented  dipoles  or  chaff  (as  in  Section  7.5) 
are  given  in  Table  2. 

Substituting  the  values  of  {SIRn}®=1  into  (7-53)  yields  an  output 
SIR  of  0.820,  which  is  less  than  the  SIR  for  one  filter  alone  (SIR  =  0.973). 
As  predicted  by  (7-55).  SIR  is  not  improved  by  using  additional  orthogonal 
filters . 


TABLE  2 


COMPUTER  RESULTS  FOR  DISTRIBUTED  PLANAR  REFLECTOR  TARGET  MODEL  IN 
RANDOMLY  ORIENTED  DIPOLE  CLUTTER, USING  SIX  ORTHOGONAL  FILTERS 


Magnitude,  I 
Phase  j 

1.516,0  1 

!  .516.  0  | 

I  ; 

1.516.0  | 

1.516,0  i 

i  i 

1.516,0  1 


Magnitude,  j  Magnitude,  i Magnitude,  'Magnitude,  Magnitude, 
Phase  Phase  j  Phase  I  Phase  Phase 


;  Phase 

i 

1  same 
las 


Phase 

same 

I  as 
!  n-l 


Magnitude,  |  Magnitude,  Magnitude,  : Magnitude.  I  Magnitude,  {Magnitude, 


.555,0,205 
.071,-2.551 
.2*5.-2.940 
.217,-2.7*7 
.020,.  1,026 


.05*, ,907 
.576.-5.155 
[  .119,-5,075 

!  .492.- .00* 

1 

.167,020 


,1*6,-2.960 

.051.2.496 

.251,-2.575 

.264,-5.0*4 

.*41,-0.029 


.075,-5.066;  .502,-0.1*2 

•356,-0.05jj  .41*,- 3.016 

.245,2  702 !  .652.0.005 

.205,2.76*  .597.5.0*2 

.248,0.410  !  .1*0,-0.145 


Horit, 

j 

Magnitude,  1 

— 

Magnitude, 

Magnitude* 

Magnitude, 

Magnitude, 

Magnitude, 

Filter 

Phase  j 

Phase 

Ph«e 

Phase 

Phase 

Phase 

f20 

.516,  0 

,535,0.015 

.066,0. 1 33 

.233,2.854 

.392,-0.254 

.131,1.696 

f2l 

.316,0 

,050.-0.595 

,546,-3.125 

.0*6,-0.222 

.355,2.995 

.138,2.819 

f22 

516, 0  j 

.491,2.86* 

.165,-3,074 

.166,2.361 

.454,-0.7*3 

.117,-0.691 

*25 

.316,  a 

1 

.  159.2.834 

.486.0.020 

.045,1.354 

.116,-2.017 

.040,-1,751 

f2* 

.316,0  [ 

j 

.076,-1.229 

.066,3.016 

.175,-0.625 

.454,2.635 

.259.2.9*0 
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7.7.3  Further  Results  for  a  Target  Consisting  of  Distributed  Planar 

Reflectors  in  Sea  Clutter 

The  use  of  multiple  orthogonal  filters  has  also  been  investigated  for 
detection  of  a  distributed  target  consisting  of  planar  reflectors  in  rough 
sea  clutter.  The  results  are  given  in  Table  3.  The  signal  and  the  first 
filter  in  Table  3  are  the  same  as  those  given  by  (7-26)  and  (7-27),  and 
SIR^  is  the  same  as  in  (7-28). 

Substituting  the  values  of  |  SIRn  j  into  (7-53)  yields  the 
following  overall  SIR  values  as  a  function  of  N,  the  number  of  filters 


used.  For  a  single  filter,  N=l, 

N=l;  SIR  =  SIRj  =  75.43  .  (7-61) 

For  two  filters, 

N=2;  SIR  =  72.95  .  (7-62) 

For  three  filters, 

N=3;  SIR  =  60.08  .  (7-63) 

For  four  filters, 

N-4;  SIR  =  48.46  .  (7-64) 

For  five  filters, 

N=5;  SIR  =  44.11  .  (7-65) 


Again,  SIR  decreases  even  though  better  detectability  is 
theoretically  obtained  from  (7-57)  -  (7-60)  as  more  filters  are  added. 
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table  3 


COMPUTER  RESULTS  FOR  DISTRIBUTED  PLANAR  REFLECTOR  TARGET  MODEL  IN 
ROUGH  SEA  CLUTTER,  USING  FIVE  ORTHOGONAL  FILTERS 


SIRn  and 

signal/filter 

coefficients 

n»l 

■ 

n«2 

n»3 

n«4 

n-5 

SIR 

n 

75.43 

70.48 

13.38 

1.85 

0.48 

Vert, 

Magnitude, 

Magnitude, 

Magnitude, 

Magnitude, 

Signal  i 

Phase 

Phase 

Phase 

Phase 

“10 

.021.3.142 

same 

same 

•« 

.004.3.142 

as 

as 

as 

■ 

U12 

.002,3.142 

n*i 

Q«t 

n»i 

n*l 

u13 

.004,3.142 

u14 

.009,3.142 

Horit. 

Magnitude, 

Magnitude, 

Magnitude, 

Magnitude, 

Magnitude, 

Signal 

Phue 

Phase 

Phase 

Phase 

Phase 

u20 

.212,  0 

tame 

tame 

tame 

U21 

.7*0,  0 

as 

as 

as 

u 

u22 

.371,  0 

n-J 

0*1 

n-l 

0*1 

“23 

.105.0 

u24 

•017,0 

Vert. 

Magnitude, 

Magnitude, 

Magnitude, 

Magnitude, 

Magnitude, 

Filter 

Phase 

Phase 

Phase 

Phase 

Phase 

fIO 

.006.0 

.026,2.139 

.011,3.076 

.234.1.461 

.078,-2.432 

fll 

.008.0 

.033,0.097 

.010.3.134 

.130.0.491 

.049,2.863 

f12 

.004,  0 

.026.3.014 

.009,3.112 

.313,0.762 

.949.0.000 

f13 

.003,0 

.021,0.126 

.002.3.076 

,147,2.914 

.049.-0.906 

f14 

.023.  0 

.007.1.226 

.041.3.129 

.401,2.907 

.132.-0.999 

Horit. 

Magnitude, 

Magnitude, 

Magnitude, 

Magnitude, 

Magnitude, 

FUter 

Phase 

Phase 

Phase 

Phase 

Phase 

f20 

.340.0 

.791,0.010 

.303,-3.133 

.044..025 

.012,1.123 

f2l 

.856.  0 

,040.-3.074 

.315.0.000 

.000,0.000 

.000.0.000 

1  f22 

.388,  0 

.604,-3.139 

.693,3.137 

.025,-2.535 

.Ot7, -2.989 

r» 

.001. 0 

3)63,2.987 

.007,2.990 

.714.-0.086 

.237.2.298 

!  f2< 

.000,0 

.009,1.341 

.001.0.001 

.361.-0.372 

.120.2.006 
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8.0  CONCLUSION 


Signal-to-interference  ratio  (SIR)  maximization  has  been  used  to 
obtain  an  optimum  signal-filter  pair  for  a  polarimetric  radar  when  targets 
and/or  clutter  exhibit  random  polarization  modulation.  The  results  can 
easily  be  extended  to  include  the  design  of  a  likelihood  ratio  receiver  for  the 
same  problem.  Considerable  insight  into  the  theoretical  solutions  has  been 
obtained  by  implementation  and  test  of  a  computer  program  to  yield  the 
maximum  SIR  and  Bayesian  systems,  i.e.,  the  "best"  signal  and  receiver 
configurations  in  each  case. 

Relatively  simple  expressions  for  the  polarimetric  scattering  function 
of  randomly  oriented  dipoles  have  yielded  expressions  for  SIR  in  some 
simple  but  important  cases,  and  these  expressions  have  been  analyzed  in 
order  to  interpret  the  computational  results.  Some  important  insights 
have  been  obtained  from  the  SIR  expression  for  distributed  planar  targets 
and  randomly  oriented  dipole  clutter,  i.e.,  for  the  typical  "target  in 
chaff"  problem.  These  insights  have  resulted  in  the  design  of  a  new 
polarimetric  clutter  canceller  (Figure  7-7  and  Eqa,  (7-36)  and  (7-42)1 
which  theoretically  allows  a  polarimetric  radar  to  "see"  through  chaff. 

The  likelihood  ratio  receiver  is  obtained  by  adding  a  set  of 
mutually  orthogonal  filters  to  the  maximum  SIR  receiver.  The  filter  output 
powers  arc  weighted  and  summed  to  implement  a  Bayes  optimum  detector 
for  a  Gaussian  signal  in  Gaussian  interference.  A  surprising  result  is 
that  each  additional  filter  reduces  overall  SIR,  even  though  detection 
probability  is  increased  more  than  false  alarm  probability  Pp  for  a 
high  threshold  setting.  This  result  follows  from  the  fact  that  PD  and  Pp 
for  a  high  detection  threshold  are  dependent  upon  the  tails  of  the 
distributions  describing  the  receiver  output  for  signal  plus  interference 
and  for  interference  alone.  SIR,  on  the  other  hand,  is  the  mean  value  of 
the  output  distribution  for  signal  alone  divided  by  the  mean  value  of  the 
output  distribution  for  interference  alone.  SIR  is  thus  not  very  dependent 
upon  tail  behavior. 
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Despite  the  shortcomings  of  SIR  as  a  measure  of  receiver  performance, 
SIR  has  proven  extremely  useful  as  a  criterion  for  receiver  design.  In  fact, 
it  is  difficult  if  not  impossible  to  specify  an  optimum  receiver  and  an  associated 
optimum  radar  signal  for  the  polarimetric  case  by  direct  solution  of  a  likelihood 
ratio  formulation.  A  solution  of  the  pi'oblem  using  SIR  maximization  techniques, 
however,  has  been  demonstrated.  A  computer  program  specifies  the  signal 
and  receiver  for  maximum  SIR  and  the  Bayes  optimum  signal  and  receiver, 
when  the  polarimetric  scattering  functions  of  target  and  interference  have 
been  specified.  The  relation  of  the  polarimetric  scattering  function  to  tapped 
delay  line  filter  and  dipole  scattering  models  has  been  discussed  in  Sections 
3.1  and  4.4. 
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appendix  a 


4 


EVALUATION  OF  INTEGRALS  FOR  COMPUTING  SIR  WITH 
SIGNAL  AND  FILTER  DESCRIBED  BY  COMPLEX  FOURIER  SERIES 

AND  SINC  FUNCTIONS 


A.  1  Fourier  Series  Description 

We  would  like  to  obtain  easily  computed  expressions  for  the 
integrals  1^  in  Eq.  (5-17)  and  ^  in  Eq.  (5-18). 

The  first  integral  is 


4^+6  4>  /  2 


I,  = 
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2 

sin  ( Tnj)T) 


4>  ~  A  /  2 
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a.  2 


_ _1  “  COSJC _ _ 
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X,  =  2itT(  4>  d>  /  2 1 

r c  3 

1  1  -  cos  x 

f  (2it(k-i)  +  x|  1 2tt(p-q)  +  x) 

xl  «  2nT| 4>s-A4>/21 


dx  . 


(A3) 


Consider  two  cases.  For  Case  1,  k-2,  =  p-q.  For  Case  2,  k-l  *  p-q.  For 
Case  1.  we  have 


I 


1  -  cos  X 

(hm  *  x)2 


dx  ’ 


(A4) 


where  =  2tr( k- 2.)  =  2n(p-q).  Letting  y  =  x  +  and  using  the 
fact  that  is  an  integer  multiple  of  2 rr,  we  have 


=  L-  f 
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xi*Hkt 
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— 2 — dy 
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Integration  by  parts  then  yields 
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where 


Si(x) 


sinc(y)dy  »  ’Si(-x) 


(A7) 


for  Case  2  (k-2.  *  p-q),  we  can  use  a  Heaviside  expansion  to 
phow  that 
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Substituting  Eq.  (A8)  into  Eq.  (A3),  we  have 
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where 


Cln(x) 
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and  y  is  Euler's  constant.  The  functions  Ci(x)  and  Si(x)  can  be  computed 
from  properties  given  on  pp.  231-233  of  Abramowitz  and  Stagun's  Handbook 
of  Mathematical  Functions.  In  Eqs.  { A6)  and  (A9), 
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where 


x 2  =  2ttT [  <is  +  A  <t>/2! 


The  integral  in  Eq.  (5-18)  is 
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A-4 


Sine  Function  Description 


A. 2 


We  would  like  to  obtain  easily  computed  expressions  for  the 
integrals  Ij  in  Eq.  (5-19)  and  I2  in  Eq.  (5-20) 

The  first  integral  is 
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where  H.  ,  and  H  are  defined  as  in  Eqs.  (All)  and  (A12). 
k*.  pq 

Using  the  identity 
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and  letting  x  =  2ttBt,  we  have 
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In  Eqs.  (A19)  and  (A20) 
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Consider  two  eases.  For  Case  l,  H,  .  «  H  i.e.,  k-£  s  p-q, 
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For  Case  2,  **  HpcJ.  For  Case  1,  Eq.  (A16)  becomes 
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From  Eqs.  (A5)  and  (A6) , 
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For  Case  2,  we  can  use  the  identity  in  Eq.  (A8)  to  obtain 
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To  evaluate  Ip,  we  can  use  the  additional  identity 
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Substituting  Eqs,  (A8)  and  (A26)  into  Eq.  (A20),  we  have 
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Letting  y  *  Hki  +  x  and  using  the  fact  that  equals  an  integer  multiple 
of  2*t,  the  integrals  in  Eq.  (A27)  can  be  written 
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It  follows  that 
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Equation  (A24)  for  Case  1,  together  with  Eqs.  (A25)  and  (A30)  for 
Case  2,  yield  the  desired  expressions  for  the  integral  1^  in  Eq.  (5-19). 
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APPENDIX  B 


DESCRIPTION  OF  COMPUTER  INPUT  /OUTPUT 
FOR  THE  SIMPLEX  SIR  MAXIMIZATION  PROGRAM 

* 

The  input  to  the  computer  program  is  given  by  E [ b-.  ^mn  1 ^  1 » 

Etc.,  c*  |t, $)] ,  p.  (4>,x),  and  p  ($,x).  These  quantities  determine  the 
ij  mn  d  c 

target  and  clutter  polarimeiric  scattering  functions  as  in  Eq.  (3-12),  where 
the  indices  i,  j,  m,  n  are  each  equal  to  one  (vertical)  or  two  (horizontal). 

The  output  of  the  program  is  divided  into  a  separate  section  for 
each  type  of  parameter  to  be  optimized.  In  Table  Bl,  the  variable  parameters 
are  the  magnitudes  of  the  signal  frequency  components.  The  table  shows 
the  eleven  vertices  of  the  initial  simplex.  At  each  vertex,  a  different  set 
of  signal  magnitudes  is  used.  In  the  example,  all  or  nearly  all  the  frequency 
components  (#0.  #1,  #2.  #3,  and  #4)  have  the  same  magnitudes,  but  this 
need  not  be  true  in  general.  It  is  only  necessary  that  the  initial  simplex 
vertices  span  the  space  of  all  desired  solutions.  Since  only  signal  magnitude 
is  varied,  all  other  parameters  in  the  table  are  fixed. 

The  best  magnitudes  for  the  five  signal  frequency  components  are 
found  with  all  other  parameters  held  fixed.  These  magnitude  values  are 
entered  at  the  bottom  of  Table  Bi,  under  the  heading  “Final  Component 
Coefficients  for  Iteration  *1.M 

The  next  output  section  (Tabic  B2)  has  the  heading:  "Maximizing 
SIR  by  Varying  Filter  Magnitude."  The  variable  parameters  are  now  the 
magnitudes  of  the  filter  components.  Since  only  filter  magnitude  is  varied, 
all  other  parameters  in  Table  B2  are  fixed.  The  best  solution  from  Table  Bl 
is  carried  down  into  Table  B2  as  the  initial  vertex  (Point  #1).  This  carry- 
down  procedure  is  used  throughout  the  program,  allowing  the  algorithm  to 
"build"  on  optimized  parameters  and  ensuring  that  SIR  is  nondecreasing. 


The  best  filter  magnitude  parameters  are  given  under  the  heading,  "Final 
Component  Coefficients  for  Iteration  #1. 

The  iteration  number  has  not  changed  because  an  "iteration"  is 
defined  as  one  cycle  through  all  the  different  parameter  types:  signal 
magnitude,  filter  magnitude,  signal  phase,  and  filter  phase.  The  last 
parameter  variation  for  Iteration  #1,  "Maximizing  SIR  by  Varying  Filter 
Phase,"  is  shown  in  Table  B3. 

The  final  coefficients  at  the  bottom  of  Table  B  3  comprise  the 
starting  set  for  Iteration  #2,  which  again  adjusts  all  parameter  types  in 
turn.  Table  B4  shows  the  final  parameter  variation  (filter  phase)  for 
Iteration  #2.  Table B 5  shows  the  same  part  of  Iteration  #3.  When  the  SIR 
improvement  between  succeeding  iterations  is  less  than  a  small  number  (0,1 
in  this  case),  the  optimization  program  is  terminated.  The  coefficients  and 
associated  SIR  at  the  bottom  of  Table  B5  then  represent  the  final  result. 

The  algorithm  may  be  sensitive  to  initial  conditions.  One  way  to 
test  for  such  sensitivity  is  to  optimize  coefficients  in  a  different  order,  e.g.  * 
by  starting  with  signal  phase  rather  than  signal  magnitude.  The  first  and 
last  parameter  variations  obtained  by  optimizing  the  phase  first  are  shown 
in  Tables  B6  and  B7.  The  results  in  this  case  are  slightly  different,  but 
the  SIP.  values  are  nearly  identical  (compare  Tables  BS  and  B7). 
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TABLE  B2 

MAXIMIZING  SIR  BY  VARYING  FILTER  MAGNITUDE 
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0  1*1X7*  i  8  137 

• 

1 

o  33*si*i 

3 

143 

9  I3H«*1  9  CM 

0  333X1*1 

0 

090 

0  I£W7*1  J  137 

• 

4 

0  f*ll*4j 

3 

142 

0  l*IH*i 

9  030 

9 

14117*1 

9 

090 

0  0*7(7* 1  a  137 

1  *73 

t 

0 

0  3|7«*4j 

3 

143 

0  1*30*1 

0  030 

0 

:o7»i»i 

0 

O'K) 

0  II«X»J  3  334 

4 

t 

0  *1  lt«»l 

3 

140 

9  I3**4*i 

0  930 

9 

314X74, 

0 

090 

0  333X2*1  3  £34 

3 

3 

o  i**e*»i 

3 

143 

0  3*30*1 

0  C30 

0 

339X7*1 

0 

coo 

0  l*lf**j  3  334 

4 

3 

0  39 ««*»i 

3 

143 

0  13314*1 

9  039 

0 

333X7* , 

0 

030 

0  i:c*2*i  3  334 

4 

4 

0  SBlXi*, 

3 

142 

0  l*lU*| 

0  030 

0 

>41*7*1 

9 

090 

0  0*7tl*i  3  £34 

1  <33 

to 

0 

0  31  7ft* | 

3 

149 

9  |43XI*| 

0  039 

0 

307X7*1 

9 

000 

O  1IHl*l  4  314 

to 

1 

0  4| |(>4 i 

3 

143 

0  13*0*1 

0  939 

0  3ta(3*i 

0 

oco 

0  333(7*1  4  3l« 

to 

a 

0  &»*Xt#| 

3 

142 

9  1*30*1 

9  930 

0 

333X3*1 

0 

000 

0  I*IX«*1  4  314 

to 

3 

0  33*04  j 

3 

142 

0  1390*1 

0  030 

0 

333X1*1 

0 

000 

0  l£3f 7* i  4  314 

10 

4 

0  70104  4 

3 

142 

0  141X4*1 

9  039 

0 

mty*i 

0 

000 

0  0*7X1* 1  *  314 

3  033 

It 

0 

0  31*044 

3 

142 

9  I43XX* i 

0  039 

0 

397X7*  j 

0  009 

0  1 1*47* 1  3  3*3 

It 

l 

0  41  104  j 

3 

142 

0  13*0*1 

9  030 

0 

314X3*1 

0,  009 

0  333X7*1  3  3*3 

It 

a 

0  9**04  4 

3 

142 

0  £»30», 

0  030 

0 

333XI*i 

0 

099 

0  t*l(l»j  J  3*3 

tl 

3 

0  31*04  j 

3 

142 

0  1 310*  , 

9  030 

0 

333X(*| 

0  909 

0  1 £4£l*i  9  3*3 

It 

4 

0  £0104  4 

3 

142 

0  I4|«*, 

0  039 

0 

141X1*1 

0 

990 

0  0*7£7*i  3.  3*3 

3  4*1 

*tnal  C»**ana*t  c»»**uit"««  *<r  Ilium*  •  I  . 


anal  Ca»a  — — - - . filna 


• 

V«rt 

ttl 

M«r  1 1 

#fl 

van 

*al 

Marti 

*•1. 

Oft 

0 

0 

31  Cttj 

3 

i*a 

0 

t43fl») 

0  090 

0 

307X7*1 

0 

oco 

0 

ll*X»l>! 

0  000 

Of  8 

l 

0 

41  l«it j 

2 

|4J 

0 

IJtft#, 

0  coo 

0 

1iafX*j 

0 

oco 

0 

33X7*1 

0  OOO 

On 

a 

0 

Sttflf , 

a 

1*3 

0 

ataiit , 

0  coo 

0 

MJEJl*  , 

0 

oco 

0 

1*147*1 

0  000 

Oft 

3 

0 

3*t<fP  4 

2 

1*3 

0 

0  coo 

0 

333(14  t 

0 

ooo 

0 

134X1*  1 

0  ooo 

Oft 

4 

0 

24MCI 

9 

1*3 

0 

l4lUt} 

0  090 

0 

141CJ*  | 

0 

000 

0 

0*  7(4*1 

0.000 

31* 


a  »« 


TABLE  B5 


MAXIMIZING  SIR  BY  VARYING  ONLY  FILTER  PHASE 


3tartln|  Caapanan*  CiiMtOinlt  *or  Itarattan  *  3  '■ 


Paint  Caa*. 
a  Prana 


1 

1 

1 

l 

1 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 

3 

3 

3 

9 

3 

a 

4 
» 
4 
4 


7 

7 

• 

• 

a 

a 

a 

a 

a 

a 

■:o 

10 

to 

to 

10 

II 

11 

II 

II 

It 


0.  31  7EX8  j  3.  143 
0  41  1EX8)  3.  143 
o  3**ex8)  3.  143 
0  334EXP j  3.  143 
o.  aatcxPj  3.  its 

0  3lTtXPj  3.  143 
0  41 lEXPj  3  143 
0.  544EXP  j  3.  143 
0.  334EXP  j  3.  143 
0.  SaiEXPj  3.  143 


0  0.  317EXPj  3.  143 
1  0.  41  lEXPj  3.  143 

3  0.  344EXPJ  3.  143 

3  0.  394EXPJ  3,  1*3 

4  0  3» lEXPj  3.  143 

0  0.  317EXPJ  3.  143 
1  0*11818)3.143 

a  o  jaacxpj  3.  M3 

3  0.  394€XPj  3  143 

4  0.  3* lEXPj  3.  143 

0  0.  317CXPJ  3.  143 

1  0.  .1  lEXPj  3.  143 

3  0  SaaEXPj  3  143 

3  0  334EXPJ  3.  143 

4  o.  aaifxpj  3. 14a 

0  0.317818)3.143 
1  C  It lEXPj  3  143 
3  0  3**818)  3  143 

3  0  3J4C»Pj  3.  143 

4  o.  aaitxp,  3. 143 

0  0.  3l7fXPj  3.  143 
1  y.  41  |(IP  j  3.  143 

a  o.  oa«xPj  3  i4u 

3  o  ajaaxpj  3. i*s 

4  0.  atisxPj  3. 143 

0  0.3171X3,3.14' 
I  0.  411(18’  9. 1*3 

3  0  JI«W|  3.  143 

5  o  aaacxPj  a.  us 

4  0.  3313X3  j  3.  143 

Q  0  il7«I»j  3.  t '» 
1  0  41  lUPj  3.148 

i  0.  jaafiPj  3. 143 

3  0.  39*8*8)  3.  143 

4  0,  34  It  IP  j  3.  143 

0  0  317818)  1  i*S 
I  0*1 'EXP j  3  l«3 

3  3.  9**8 *8)  3  143 
a  0  33*818.  3.  143 

4  0  S»UxPj  3.  1*3 

0  0  317818)  3  1*3 

I  0  4IKXPj  3  143 

3  0  3**4'*8)  3  1*3 

3  0  33*8*8)  3  1*1 

4  0  34  KIP  j  3.  143 


0  143EXP)  0  030 
0  139EXP  j  0  000 
0.  393EXP  1  0.  030 
0.  138EXP  j  0.  030 
0.  14KXP  1  0.  000 

0  143CXPj  0  000 
0. 134E«P j  0  030 
0.  343EXP  j  0.  000 
0  138EXP  j  0.000 
0.  141C)fPj  0.  COO 

0.  143E.JP  1  0  030 
0.  !3*VXPj  0.  000 
0.  343EXP)  0.  COO 
0  138EXP  j  0.030 
0.  14tCXPj  0.  OOO 

0  143EXP  J  0  030 
0  139(ltPj  0.030 

0.  aaotxp  j  0  000 
0.  IOBEXPj  0  030 
0.  1418*8)  0.000 

0.  1*38X8)  0.  000 
0  I04CXP j  0.  000 
0.  343EXP  j  0.  000 
0  IOBEXPj  0.030 
0  I41CXPJ  0.000 

0  143IXPJ  0.000 

o  toattxPi  o  ooo 
o  aa-jfXP  J  o.  ooo 
0  134EXPJ  0.030 
0.  1418*8)  0.000 

0  1 4  IMP  |  0  OOO 

o.  taatJpj  o.ooo 
o.  aaaitPj  o  coo 

0  1388X8  J  0.000 
0.  1*18X8)  0.090 

0.1*38X8)  0,030 
0.  i 3*8*8 j  4.000 
0  34341*  J  0.  OOO 
0.  ‘38818)  0  OOO 
0.  141818)  0.000 

u.  i*3cxpi  o  ooe 
0  ijatxPi  o  ooo 
o. aaj  tpj  o.ooo 

0. I34BXP  J  030 

o.  i*iF»Pj  e.000 

0  laXIP  i  0  OOO 
0  toatppj  o  03 <t 
y  a*3f»P,  o  ooo 
0  1 341 IP  ’  0  030 
0  14KXP >  0  030 

3  ItOClPj  0  ooo 
0  13*8"  0  000 

0  3*3(X*>  0  030 
0  I 38EXP i  0  OOO 

o.  lutxP  t  o  ooo 


0,  007LXP  j  0  ooo 
0  314EXP)  0,  OOO 
0  533EXP  J  0  OOO 
0  050EJCP  J  0  000 
0.  14KXP  j  0  000 

0  307EXP  j  l  047 
0.  3148*8)  l  047 
0.  3SSEX8)  1.  047 
0.  033EXP  j  1.  0*7 
0.  141EXP  J  I.  047 

0. 307EXPj  3  044 
0.  SttEXPj  3.  044 
0  SSJEXPj  3.044 
0.  300CXP  j  3. 044 
o.  ukxpj  a.  oaa 

0.  307CXPJ  3.  143 
0.  3I4EXP  J  3.  1*3 
0.  OOSEXPj  3.  143 
0.  330EXP  j  3.  143 
0.  14ICXPJ  3.  143 

0  30703*  J  4.  114 
0.  SltEXPj  4.  184 
0.  SSXlP)  4.  184 
0.  390EXP  j  4.  184 
0.  14KXP j  4.  184 

0.  307**8 J  9  334 
0.  JUEXPj  3.  334 
0.  93X18)  S  334 
0.3938*8)  9.334 
0.  14KXP j  9  334 

0.3078XPJ  0.000 
0  314X^1  0.  COO 
0.  93X18)  0.  OCO 
0.  3333  XP  j  0.000 
0.  1418X8)  0.  OOO 

0.  3078X8  1  0  OOO 
0.  9148*8)  0  OOO 
0  9938X8  )  0.000 
0.  333818  )  0  OCO 
0.  1418X8)  0  OOO 

0  307f *8)  0  OOO 
0.  5148X8 )  0.  COO 
0  1938X8  )  0  OOO 
0,  333818  )  0  OOO 
0.  1418X8)  0.  (WO 

o  3078X8)  0  050 
0  314418  )  0  OOO 
0  3338XP)  0  CM 
0  3338X4 1  0  000 
0  I4KX8)  9  030 

0  OOTfpP)  0  OOO 
0  314818  )  0  OOO 
0  3338*8 (  9  COO 
0  3311*8 )  0  OOO 
0  14KIP)  0  COO 


818 


3.447 


o  tiaexP)  o  ooo 
0.3338*8)  0.000 

o  taiexPj  o  ooo 

0  tSSEXP j  0.000 
0  04  7CXP j  0.  OOO 


o  ltacxPj  o  ooo 

0  335CXP  j  0.  000 
0  14 lEXPj  0  OOO 
0  128EXP j  0.000 
0  0478X8)  0  OOO  3.413 

0  I14CXP)  0.000 
0.  3338X8)  o  OOO 
0  14KXP  j  0  OOO 
0  1388X8)  0.000 
0  0*7818)  0.000  1.443 

0.  ItafXPj  0.000 

0  3338X8  )  0.  OOO 

o.  mcxPj  o.  ooo 

0  ISCCIPj  0.000 
0  0*7818  j  0.  OOO  1.  834 

o  iiac*Pj  o  ooo 

0  333818)  0  OOO 
0  14I8IP)  0  OOO 
0  1388 IP  J  0.  OOO 
0.0478X8)  0.000  1.443 

0.  114SXP)  0.000 
0  3338X8  )  0  OOO 
o  14I8XP)  o  ooo 
0  1388X8)  0  000 
0  047CXP)  0.000  3.413 

0  1148X8)  1,074 
0  333818)  1.074 
0  taKXPj  1,074 
0  1388X8)  1.074 
0  0478X8)  1.074-  3.344 

0  1148*8)  3.  137 
03338*8)  3.137 
0  14)8*8)  3.  1371 
0.  1388*8  )  3.  197' 

0.  047*18  j  3.1971  1.479 

0  1148*8  )  3.834 
0  33X18  )  3.334 
0  141818)  3.934 
0  1388*8)  3334, 

&  0478X8)  3.334,  1.833 

0  1148X8)  4  .114 
0  33X18 )  *  314 
0  !4K*8)  4  314 
0  1388*8)  4  914 
0  0478*8)  4.  314  3.033 

A  u*8l8j  3  343 
0  3338*8)  9  343 
0  1 4K*8  j  9  943 
9  138X18)  3  343 
0  047888  )  3  343  3  481 


81**1  (*•****«!  Ca«**itU*tt  **!•  !*****••*  •  3 


8*t*4 

Caaa 

•Nti4aa**i 

*•** 

.  »  mmm 

•■'*  *  ** 

--*  F  l  1  Iff—’ - 

.  ft 

8»**4 

V#r<4 

fftl. 

Marti 

8*1 

v#m 

ftftl 

Mftrn 

ftftl. 

f?M 

0 

0  31715X8) 

3 

0 

143(18 | 

0  ooo 

0  JC'lXPi 

0  OOO 

0  11*1*8) 

0 

ooo 

Oil. 

1 

0  41 18*8) 

n 

1*3 

0 

134(18) 

0  ooo 

0  9t*I<8) 

o  o:o 

0  0308*8) 

0 

ooo 

3 

0  9MCtP| 

3 

H2 

0 

3*31*8 | 

0  030 

0  !I!EE7| 

0  coo 

0  l*lf/»j 

0 

OUJ 

Oil 

3 

0  33*8*8  | 

3 

1«3 

0 

nui"i 

0  coo 

0  33X18) 

0  050 

0  1  *8818) 

•J 

ooo 

«♦» 

4 

0.  34 18*8  j 

3 

143 

0 

l«;tX8) 

0  000 

0  141(18) 

0  OCO 

0  0*rt*8 ) 

0 

ooo 

818 


3  447 


B-7 


TABLE  B6 


MAXIMIZING  SIR  BY  VARYING  ONLY  SIGNAL  PHASE 


SUrtiitf  c*«**n«nt  Coifftcirntt  for  Ilorotlon  •  1  •' 


Point  Co»». 

•  fritl 


- - FI  It  or- 

Vort.  Pol 


-  ltd 

Horn.  Pol. 


1 

l 

l 

1 

l 

a 

a 

a 

a 

a 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 

3 

3 

3 

3 

3 

4 
4 
4 
4 
* 

7 

7 

7 

7 

7 

■ 

t 

■ 

* 

4 

4 

4 

4 


0.  31AE17)  0.000 
0.  314EIPJ  0  000 
0  OloESPj  0  000 
0.  31407)  0.  000 
O.  31*07)  0.000 


0  O  314*17)  l  047 

1  0  31*07)  l.  047 

2  0  31407)  l  047 

1  0.31407)  1047 

4  0.314**7)  !  047 

0  0.31407)3.044 

1  0.31407)  2.044 

2  0.  314EIPJ  2  044 

3  0.31407)  2.094 

4  0.31407)  a.  044 

0  0  3t4fXPj  3  142 
t  0.31407)  3.142 

2  0.  31  4EXP  j  3.14a 

3  0.31407)  3-  t*a 

4  0.  31407)  3.  14a 

0  0.  3I4EXPJ  4.  1(9 

t  0.  314CXPJ  4.  1(4 
a  0.  31  4CXP  j  4.  114 

3  0.314**7)  4.1(4 

4  0.  314CXPJ  4.  1(9 

0  0.31407)  3  234 
1  0  314**7)  3  234 

a  0. 31407)  3.234 

3  0.  31 4CXP  j  3.234 

4  0.  314*17)  3  234 

0  0.31*07)  0  000 

1  0. 314**7)  0.000 

2  0.  314EXP)  0.000 

3  0.  314EXP;  0.000 

4  0.  31  4*  IP  J  0  .000 

0  0.  314(17)  0  000 
1  0.  3t***»)  0.000 
I  0  31407)  0.000 

3  0.31407)  0.000 

4  0.  3141*7  j  0.040 

0  0.  314(17,  e.  900 
1  0.  3l***»)  0  rtO 
3  0.  914**7)  0  000 

3  0.  914(17)  0,000 

4  0.  314017 1  0.000 


0  314EXP|  0  030 
0  31SE17 i  0  OSO 
0  316**7  J  0  C30 
0  314£x»  1  0  OOO 
0  316EXP  )  0.  030 

0  3UEIP)  0  030 
0  314EXPJ  0.  030 
0.  314EXP  j  0.  030 
0  314EXP j  0  COO 
0.  314CXP  j  0  000 

0.  314EXP )  3.030 
0  314EXP j  0  030 
0  314EXPJ  0  030 
0  314EXP  j  0.  000 
0.  314**7 I  0.  030 

0.  314CXP  j  0.  030 
0.  314EXP  j  0  030 
0.  314CXP  ,  0.  030 
0,  314CXP  j  0.  030 
0.  31407  j  0.000 

0.  314CXP  j  0.  030 
0.  314CXP)  0.  030 
0  314EXP  j  0.  COO 
0.  314E1PJ  0.000 
0.  314EIP  j  C.  000 

0  314*17  |  0.000 
0.  314EXP  j  0.  000 
0  314(17)  0.000 
0.  314**7  j  0.  030 
4  31*117  )  0.  000 

0.314*17)  1.074 
0.  314(17)  1.074 
0  314*17)  1.074 
0  314*17)  1.074 
0.  314(17)  1.074 

0.  314**7)  2.  1ST 
0  314*17)  t.  137 
'  0.314**7)  I,  137 
0  314*17)  1. 117 
0.  314*17)  I.  117 

0.  314*17)  3  234 
0  314*17)  9  194 
0.314(17)  9.194 
0  31UR7)  9.934 
0.  314*47  J  9. 134 


0.  316**7)  0  000 
0.  3t4E*P)  0  000 
0  316EXP)  0  000 
0  3UEXP)  0  000 
0.  314EXP)  0.  000 

0  3t*EXP)  0  000 
0  314SXP)  0.  000 
0  31607,  0.  000 
0.  314CXP)  0  000 
0.  314CXP)  0.  000 

0  314CXP )  0  000 
0.  314(17  )  0  000 
0  314EXP)  0.  000 
0.  314**7)  0.  000 
0.  314CX7)  0.  000 

0.  314(17)  0.  000 
0  314**7 )  0  000 
0.  314(17)  0  000 
0.  314(77  )  0.  000 
0.  314(17)  0  000 

0.  314(17  )  0  000 
0.314(17)  0  000 
0.  314(17)  0.  000 
0.  314(17)  0  000 
0.  314*17)  0.  000 

0.  314**7  J  0  000 
0  31**17 |  0  000 
0.  314(17)  0  000 
0.  314(17)  0  000 
0.  314(17)  0.  000 

0.314(17  )  0  000 
0.  314*X7  )  0  000 
0.  314(17  )  0  000 
0.314(17)  0.000 
0.  3141*7)  0.000 

0.  314«I7|  0  000 
0.314(77  )  0  000 
0.  314*7*)  0.  000 
0.  314**»  J  0.000 
0  314**7)  o  000 

0„ 914*17  )  0  000 
0.  914*1*)  0  000 
0.  914*17)  0  000 
0.  314* (7)  0.  000 
0.  9144117)  0.  000 


0  314EX7)  0.  OOO 
0  314CX7)  0  000 
0  314EX7)  0  000 
0316*17)  0.000 
0.  3UEX7  )  0.000  0.  44* 

0  314(17)  0.  000 
0  316**7)  0.000 
0  314(17)  0.  000 
0.  314EX7  )  0  000 
0.314*17  )  0.000  0.734 

0.314*17)  0.000 
0  31*5X7)  0  OOO 
0.  314(17)  0.  OOO 
0.  31*(X7)  0.000 
0.  3145X7)  0.  OOO  1.13a 

0  314(17)  0.000 
0.  314(17)  0.000 
0.  314**7)  0.000 
0.  314(17)  0.000 
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APPENDIX  C 


ON  MULTI-CHANNEL  DETECTION  OF 
RANDOM  SIGNAL  IN  GAUSSIAN  NOISE 


ON  MULTI-CHANNEL  DETECTION  OF  RANDOM 
SIGNALS  IN  GAUSSIAN  NOISE 


1.  INTRODUCTION. 

In  the  classical  books  Van  Trees  [l],  [2]  has  presented  the  single  channel 
detection  theory  when  signals  are  deterministic  or  have  random  parameters. 
Extension  of  these  works  to  multi-channels  case  has  received  little  attention. 
An  earlier  reference  is  [3]  in  which  Lindsey  derived  the  optimum  receiver  and 
its  performance  for  independent  Rician  fading  multi-channels. 

Vannicola  [4]  introduced  a  natural  extension  of  the  above  single  channel 
model  for  the  two-channel  radar  detection  problems.  In  his  model,  Vanni¬ 
cola  used  the  scatter  matrix  B  to  describe  the  characteristics  of  each  channel 
and  the  correlations  between  two  channels.  These  models  are  important  for 
both  radar  detection  and  communication  problems.  Our  revisit  to  Vanni- 
cola’s  model  reveals  that,  for  some  of  the  important  detection  problems,  the 
extension  of  single  channel  results  to  multi-channel  situations  is  not  straight¬ 
forward  and  requires  additional  techniques  to  yield  desired  expression  for  the 
optimum  receiver. 

In  this  paper,  we  consider  one  of  the  important  cases  of  such  multi- 
channels  problems,  i.e.,  the  slowly  fluctuating  point  target  model  when  the 
noise  is  white  or  colored.  Using  the  maximum  likelihood  criterion,  orthonor- 
malizing  procedure,  and  eigen  value-eigen  vector  approach,  we  have  derived 
the  optimum  receiver  for  the  two-channels  case.  Extension  of  these  results 
to  more  than  two  channels  is  straightforward. 
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2.  OPTIMUM  RECEIVER  FOR  SINGLE  CHANNEL  CASE. 

This  section  contains  a  brief  summary  of  single-channel  detection  results. 
These  results  are  well  known  and  most  conveniently  available  in  Van  Trees. 
The  reason  to  reproduce  some  of  these  results  here  in  that  it  allows  an  easy 
reference  as  well  as  a  comparison  with  multi-channel  results  derived  in  the 
remaining  sections. 

For  the  detection  of  slowly  fluctuating  point  targets  in  the  presence  of 
additive  white  Gaussian  noise,  Van  Trees  derived  the  statistical  model  as 

f(t)  =  y/Ft  bf(t)  +  i5(t)  0<t<T:Hx  (2.1) 

f(t)  =  u>{t)  0<t<T:H0  (2.2) 

where  f  (f)  is  the  complex  envelope  of  the  received  wave  form,  6  is  a  zero-mean 
complex  Gaussian  random  variable  which  satisfies 

£(|6|J)  =  £(66*)  =  2  a?,  (2.3) 

f(t)  is  the  complex  envelope  of  the  transmitted  signal  with  unit  energy,  i.e., 

f  irwi’* = /  mrm = i,  (2.4) 

w(t)  is  a  zero-mean  complex  Gaussian  white  noise  process  which  is  indepen¬ 
dent  of  b  and  satisfies 

£(w(0w*(tt))  =  JW(c-u).  (2-5) 

The  optimum  receiver  is  constructed  by  using  the  complete  orthonormal 

(C.  0.  N. )  set  expansion,  and  by  finding  the  sufficient  statistic.  It  turns  out 
that  this  sufficient  statistic  has  the  form 

ri=  f  f(t)f-(t)dt.  (2.6) 

Jo 
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The  optimum  receiver  computes  |i?i|2  and  decides  acceptance  or  rejection  of 
Hi  depending  upon 


H  , 

2  >  .* 


W  ?  r 
Ho 


(2.7) 


where  r*  is  called  the  threshold.  This  structure  of  the  optimum  receiver  is 
graphically  presented  in  Figure  1  below: 


Figure  1.  Optimum  receiver  of  Single-Channel  Case. 


3.  OPTIMUM  RECEIVER  FOR  MULTI-CHANNEL  MODEL. 
WHITE  NOISE  CASE. 

For  the  multi-channel  model  derivation  of  the  optimum  receiver  is 
relatively  complicated.  For  simplicity  in  presentation  we  consider  the  two- 
channel  case  first. 

The  statistical  model  for  two-channel  case  can  be  formulated  as  follows: 


1(0  =  4/(0 

+  w(f) 

0  <  t  <  T 

:  Si 

(3.1) 

£(0  = 

2t(0 

0<  t  <  T 

:  Ho 

(3.2) 

where  the  vectors  r ,  /  and  tv  and  matrix  6  can  be  expressed  in  terms  of  their 


components  as  follows: 


1(0 


&(0 


U(0j 

,  /to  = 

' /.(O' 
L/iWj  ’ 

f  u»4(e) 

.  i  = 

bn  b \2 

|_u/2(e) 

bn 

(3.3) 
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The  vector  f(f)  is  the  convex  envelope  of  the  received  wave  forms,  b  is  the 
scatter  matrix  of  the  target,  f[t)  is  the  convex  envelope  of  the  transmitted 
signals  and  its  energy  equals: 

e=  fmfimydt  =  f iawi1^  +■  /Vw  iJ* 

Jo  ~  ~  Jo  Jo 

=  E{  +  E2  (3.4) 

where  E\  and  E2  are  the  integrals  of  |/i(0l2  and  |/a(0!2  an<*  represent 
energies  in  channel  1  and  2,  respectively.  (For  amplitude-modulated  signals 
actual  transmitted  energy  equals  E/2).  The  zero-mean  vector  white  noise 
process  xb(t)  is  independent  of  the  scatter  matrix  l  and  the  two  components 
of  the  vector  w{t)  are  independent  of  each  other.  Moreover,  it  is  assumed 
that  the  two  components  have  equal  spectral  densities  in  the  two  channels, 
i.e., 

*!*(0&*MI“Wo  \S[‘oU)  s ri-4 

■  4  1 

=  N06{t  -  tt)/.  (3.5) 


Following  the  procedure  of  deriving  the  optimum  receiver  for  the  single¬ 
channel  model,  we  define  a  vector  C.  0.  N.  set  {j^(*)}n»t*  Since  this  is  a 
C.  0.  N.  set,  the  elements  £.(t)  satisfy 


^i«(0 

^a»(0 


and 


rT 

~}  Jo 


(3.6) 
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where  t  denotes  complex  conjugate  and  matrix  transpose.  In  terms  of  ^.(t)’s 
the  received  vector  r(£)  cam  be  written  as 

k 

i{t)  =  .lim  (3-7) 

i  —  l 

where 

fi  =  [  iT(tw m  =  l  (m 

Jo  Jo  ~ 

From  the  assumptions,  made  earlier,  on  the  process  r(t)  it  is  easily  observed 
that  the  coefficients  Vs  are  zero-mean  complex  Gaussian  random  variables. 
If  we  can  choose  the  C.  O.  N.  set  such  that  only  k  of  the  coefficients  are 
dependent  on  which  hypothesis  is  true  and  if  these  k  coefficients  are  statis¬ 
tically  independent  of  the  other  coefficients,  collectively  they  will  give  us  a 
^-dimensional  sufficient  statistic.  The  coefficients  ?<'s  will  satisfy  the  inde¬ 
pendence  condition  if 

=0  for  l  =  0, 1  and  for  all  *  £  j  (3.9) 

On  Hq  the  desired  condition  is  satisfied  due  to  the  fact  that  r<*s  depend 
only  on  the  white  noise  process.  On  H[,  the  correlation  between  V s  can  be 
written  as: 

=  JJ  ^  +  w,,  (3-10) 

where 

£,(!.«)  =  £[i/(i)r  Mi' I  (3.11) 

is  the  kernel  of  received  signal  components.  This  kernel  is  known  provided 
all  statistics  of  the  scatter  matrix  6  are  known  and  the  transmitted  signal 


wave  form  f[t)  are  given.  From  (3.10),  f,’s  are  uncorrelated  if  and  only  if 
the  functions  in  the  C.  0.  N.  set  satisfy  the  integral  equation 


J I  -  A«A;. 

or  equivalently 

J  =  A^.(t).  (3.12) 

This  last  expression  shows  that  if  the  eigen  functions  4>i[t)'s  and  eigen  val¬ 
ues  A;’s  can  be  solved,  then  the  construction  of  the  likelihood  ratio  test  is 
straightforward.  In  other  words,  the  likelihood  ratio  is  given  by 


A(J2)  =  lim 

v  1  k—oc 


P{Rk\Hi\ 


w. 


Urn 


(3.13) 


The  probability  density  functions  of  the  complex  Gaussian  variables  as 
appears  in  (3.13)  are  given  by 


=  MM,  +  Xd)-‘«*i>{-|A|V(Mo  +  *<)) 


(3.14) 


Substitution  of  (3.14)  in  (3.13)  and  then  taking  its  logarithm  and  further 
simplification  results  in  the  following  log  likelihood  ratio  test 

'  A, 


In  —  lim  V  —p-—\Rt\2  >  i 


(3.15) 


where  -y  «  In  q  4*  £  In  {(ty,  +  A,)//V0}  is  the  decision  threshold.  If  the 
number  of  nonzero  A,'s  is  infinite,  then  construction  of  the  optimum  receiver 
according  to  (3.15)  is  impractical,  and  instead  we  must  find  its  closed  form 
representation.  For  the  single-channel  model  it  was  straightforward  to  solve 
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(3.12)  because  the  function  /( u)  can  be  separated  from  the  kernel,  i.e.,  in 
the  single-channel  case 


K.(  t,u)  =  £(B/(i)/(«)5-|  =  E(  |ip)/(t)/-(u).  (3.16) 


Due  to  this  simplification  the  integral  equation  (3.12)  can  be  written  as 

•  r 
'0 


E(\b\2)f(t)  [T  /»*,( u)du  =  A ,-*(*).  (3.17) 

Jo 


It  follows  that  a  solution  of  (3.17)  is  given  by  ^(t)  =  /(0/\/*t  ant^  the 
associated  eigen  value  =  et£(|6|J)  =  2o£ee  =  er.  Any  other  function 

0«(t)»*  =  2,3,...  which  is  orthogonal  to  <i>i (t)  is  a  solution  of  (3.17)  with 
zero  eigen  value.  Thus,  in  brief,  the  single-channel  case  is  straightforward 
due  to  aforementioned  simplification  and  the  corresponding  log  likelihood 
ratio  statistic  is 


W0(W0  +  «.) 

which  can  also  be  written  as  (2.7). 

For  the  multi-channel  case,  in  general,  the  kernel  /£({,«)  cannot  be  writ¬ 
ten  in  the  form  of  (3.16).  Hence  we  cannot  derive  the  sufficient  statistics  in 
the  same  way  as  in  the  single-channel  case.  In  this  paper,  we  first  consider 
some  special  cases  which  lend  themselves  to  easier  solutions.  Our  next  goal 
is  to  show  that  it  is  possible  to  extend  the  results  so  obtained  to  a  gen¬ 
eral  setting  and  thus  obtain  the  optimum  detector  of  Gaussian  signal  in  the 
presence  of  white  or  colored  noise. 


C-7 


CASE  I.  Identical  Signal  Envelopes. 

Suppose  that  the  two  transmitted  signals,  in  two  different  channels,  have 
identical  envelopes  and  may  differ  only  in  their  energies  and  phases,  i.e., 

".M 

[At 


2(0  = 


hit) 

hit) 


5(0  = 


s/Fi  t>9x 

e i*' 


5(t)  =  As{t)  (3.18) 


[T\fi[t)\*dt  lm  0, 
Jo 


1 


where  E\  and  E-i  are  the  energies  and  6\  and  are  the  phases  of  hit)  and 
h (Os  respectively,  and  s(f)  is  a  real  valved  function  with  unit  energy.  Thus, 

-T 

' o 

and  (3.19) 

-r 

sr(t)dt  =  1. 

'o 

The  signal  components  in  received  wave  forms  can  be  written  as 


f 


k~ft,\ ...  [  4u  bn 
-  !  in  bn 


where 


4i(0  = 


Wi  y 

VF  t}6 

m  m 

s(0  =  xi^(0  +  *2^3  (0  (3.20) 


•it) 

0 


»  iM- 

*i  =  6n\/ E\t^1  +  »  6f  A  *  A^6i 

X2  a  621%/ E^t}Sl  +  62,  ~  62" A  ~  Ar6j 

fcl- 


611 


621 

622 


and  62  = 

We  observe  that  ^  (0  and  <£2(0  are  orthonormal  functions,  the  signal  process 
6/(t)  is  a  two-dimensional  vector  in  and  <£2  plane,  and  its  components 
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X\  and  12  are  random  variables.  These  two  functions,  d^t),  d2( 0»  can  be 
augmented  into  a  C.  0.  N.  set.  In  summary,  it  is  possible  to  expand  £(t)  as 


OO 


’  Y  -  Y  +  Y  ^,(0 :  (3-21) 


t=i 


t=i 


t=3 


£(0» 


I  oo  oo 

»=1  t=l 


:  Hn 


(3.22) 


where  it  is  important  to  note  that  for  i  >  3,  all  r,-’s  are  independent  of  the 
hypotheses  H i  and  Hq.  Because  of  the  independence  between  b  and  w(t)  the 
first  two  coefficients  n  and  rj  are  independent  of  the  remaining  coefficients 
T{  —  to;,  i  >  2.  But  ri  and  f2  are  not  necessarily  independent  of  each  other. 
In  any  case,  this  does  not  prevent  us  from  obtaining  the  likelihood  function 


A/m  _  PlRuWx] 

1  j  P[J2|ff0)  P[Rulk\HQ\' 


(3.21) 


The  joint  covariance  matrix  of  can  be  derived  as  follows: 

J5(r,rj|if0|  s  E[JJ  (u)<£t  du] 

m  Jf  $\{t)[NQ6{t  -  u)\4.(u)dt  du 

rT 

—  Nq  I  £{t)£At)dt  »  N06ij  ijm  1,2.  (3.22) 

Jo  ; 


In  a  similar  manner  and  due  to  independence  of  x<  and  u>2,  i,j  -  1,2  it  is 
observed  that 


**  E[bJ  AA'b'\  +  N0Sij.  (3.23) 
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Thus,  R  =  {RitRi)'1'  is  a  zero-mean  Gaussian  random  vector  whose  covari¬ 
ance  matrices  are  given  by 


Ka  =  .EjM'lffol  =  i V 


(3.24) 


and 

Kx  =  E[M!\Hi\  =  E[b£A'b'\  +  N0I 

(3.25) 

=  £,+  N0it 

where  j£,  —  i?[6AA^|  is  the  covariance  matrix  of  the  received  signal  com¬ 
ponents. 

The  matrix  may  ke  diagonalized  by  a  unitary  transformation  T,  i.e., 

Z'KJ*  °1  (3.26) 

1°  M 

where  T  =  (C^C,)  is  a  2  x  2  unitary  transformation,  mi  and  hi  are  eigen 
values  of  Ktt  and  C_x  and  C2  are  the  corresponding  eigen  vectors.  Using  thfc 
fundamental  theorem  of  linear  algebra  we  then  expand  K,,  K&  and  Klt  in 
terms  of  eigen  values  and  vectors  as  shown  below. 

K.  =  m£,£t  +  m£>c} 

&  =  Wo£,£l  +  N o£*fil  =  NoTIT' 

Ki  =&.  +  &>=  (m  +  Wo)£i£l  +  («  +  «b)£i£l  (3'27) 

r  1 

rr 


»  T 


Hi  4*  No  0 
0  mj  +  Wo 


The  above  results  also  allow  us  to  write  the  inverses 


££l  »T  1  K 


k  0 
0 


Tf,  =  T 


iV04-Ml 


rf.  (3.28 


W  VJ  4.  - 

Substitution  of  the  joint  Gaussian  probability  density  functions 
PIJMaW  =  (2»r) -*|*j|-1/*«a*»{-jRt  jrr‘S3  .*  =  0.1 
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in  (3.21)  and  simplification  of  its  logarithm  gives  the  two-channel  log  likeli- 

t 

hood  ratio 


ht=SH&l-KZl)R-{R'T) 


(mi+So)Wo  0 
°  lm  +  N0)No 


=  £ 


m» 


fr[  (m*  +  Nq)nq 


\Y<f 


{t'r) 

(3.29) 


where 


Y  =  T'R  =  [Cl  :  C2]t^  = 


Cj  R 
Cj  R 


(3.30) 


The  R  vector  is  derived  by  correlating  the  received  waveform  with  the  trans¬ 
mitted  signal  s(t),  i.e., 


£= (ri,r2)T  =  f  timm* 

=  [  $(t)r(t)dt.  (3.31) 

Jo 


Hence,  the  optimum  receiver  for  the  identical  envelope  case  has  the  form 
given  in  figure  2  below. 


Figure  2.  Optimum  Receiver  for  Identical  Envelope  Case 
Using  the  new  orthonormal  functions 

0j(O  3  £i3(0*  1  =  0,1 


C-ll 


gives 


rT 

Yl=  ±\{t)r{t)dt  1  =  0,1 

Jo 

and  con-jquently  the  optimum  receiver  has  the  form 


Figure  3.  Another  Form  of  Optimum  Receiver  for  Identical  Envelope  Case 


CASE  II.  The  Orthogonal  Envelopes  Case. 

This  is  another  extreme  situation  in  which  we  assume  that  the  trans¬ 
mitted  signals  in  the  two  channels  are  orthogonal.  This  implies  that 

■r 
<0 

Hence  the  signal  comonents  of  the  received  waveform  can  be  written  as 


/ am/;w*-o 

Jo 


=  6 


r*.. 

'a  (o' 

[b2 

l  ^22 

AM 

u 

7iWl 

0 

+  *ai 

0 


J  +il1  hi  +b>1  Uwl 


where 


tfi)  - 


1 

vzr 

i 


A(0 

o 

'  0 

hit) 


LAWJ 


.*,(0 


.  j_  r  Awi 
0  J  ' 
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It  is  easy  to  verify  that  <t> :[t) , <fr2{t) , ^{t) ,  and  $^(t)  form  a  set  of  orthonor¬ 
mal  functions.  They  can  be  augmented  to  a  C.O.N.  set.  Using  this  C.O.N. 
set,  it  is  possible  to  expand  r[t)  and  then  verify  that  the  coefficients  r»’ s 
are  statistically  independent  random  variables  for  i  >  5,  and  that  they  are 
independent  of  the  hypotheses  H\  or  Hq.  As  for  the  coefficients  ri,r2,r3, 
and  r*  it  can  be  seen  that  they  are  statistically  independent  of  r,-’s  i  >  5,  \ 
and  constitute  a  four-dimensional  sufficient  statistic. 

Let  R  =  (ri,r2, r3) r4)r.  Then  the  logarithmic  likelihood  ratio  can  be 
written  as 

H  i 

<*»£*(£?' -icr*  )&<i 

Ho 

where,  in  this  case,  both  and  K_x  are  4  x  4  matrices  and 

Ka  =  NoL  Ki  =  K*+  Ko  =  E{K  XT\+  NqI 
&-(bnV*L  hiy/El,  bu  y/E2^  hiVEi)- 
As  in  case  I,  the  statistic  can  be  written  as 

I  yf  J*  If  jy 

S  =  J-  <u*g0“1|(Ml  +  No)No'  (M2  +  No)No‘  (MJ  +  iVo )iV  («,  +  No)Nol~ 

=  £  (w  +  Wo)JVo|y;|! 

where 

y=t'r  =  [a&iSaSuH = w.  ^  n,  V4)r, 

Here,  as  in  the  previous  case,  Hi's  denote  the  eigen  values  and  C,’ s  denote  the 
corresponding  eigen  vectors  of  the  matrix.  The  optimum  receiver  for  this 
case  has  the  form  given  in  figure  4,  where  a*  =  M*/ {(m%  4-iV0)iV o},»  =  1, . . .  t4. 
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Figure  4.  Optimum  Receiver  for  the  Orthonormal  Envelopes 


It  is  useful  to  note  that  unlike  the  identical  envelope  case,  here  we  have 
a  four  dimensional  sufficient  statistic.  Both  cases  give  correlation  receivers. 
As  long  as  all  statistics  of  the  scatter  matrix  6  and  the  energies  of  the  trans¬ 
mitting  signals  are  known,  it  is  easy  to  derive  the  structure  of  the  optimum 
receiver  which  requires  solution  of  some  matrix  eigen  values  problem. 

Knowledge  gained  in  this  case  can  now  be  employed  to  obtain  the  opti¬ 
mum  receiver  in  the  general  situation.  This  case  is  considered  below. 

CASE  III.  The  General  Model. 

In  general  the  signal  components  in  returned  waveforms  can  be  expressed 
as 

4/M  =  »u  /l0(1)l  +*•»'  A0(t)  ff  +>»  /“(,)  • 

m  *■  «•  L  J  a 
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Whenever  the  additive  noise  is  colored,  a  possible  solution  is  obtained 
by  first  applying  a  whitening  filter  and  then  using  the  results  known  for 
the  white  noise  case.  We  follow  the  above  approach  to  solve  the  problem 
in  the  multi-channel  model  also.  More  specifically,  let  the  model  for  slowly 
fluctuating  point  target  in  the  presence  of  colored  noise  be  given  as  follows: 
r  kf(t)  +  &c(t)  +  rv[t)  -  bf(t)  +  n(t)  :  Hi 


£(0  = 


l  nc(«)  +  u/(f)  -  n(t)  :  H0 

Then,  we  can  design  a  whitening  filter  hw(t,  u)  such  that  after  passing 
through  this  filter,  the  noise  n(t)  will  become  a  white  noise  process  with 
a  height  of  spectral  density  1.  Thus,  if 


£.(0  =  J  ttte(u)<*tt 


then 


After  passing  through  the  whitening  filter  Aw(t,u),  the  received  waveform 
r[t)  becomes  r,  (t)  where 


bf.[t)  +  n.(t)  :  Hi 
n.[t)  :  Hq 


and 

/.(«)  =  J  Aw(«,tt)£(tt)du. 

Results  of  Case  III,  when  applied  to  r.(t),  given  above,  give  us  the  desired 
optimum  receiver.  The  optimum  receiver  has  exactly  the  same  form  as  given 
in  figure  4  except  that  all  expressions  are  replaced  by  their  starred  versions 
which  identifies  that  we  are  dealing  with  the  filtered  process.  The  whiten¬ 
ing  filter  Aw(t,u)  is  obtained  by  solving  a  matrix  eigen  value  problem  that 
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depends  on  the  statistics  of  the  nc(f)  process.  Hence,  the  colored  noise  case 
can  be  handled  only  if  statistics  of  the  colored  noise  are  known.  The  nature 
of  the  equation  to  obtain  hw(t,u)  remains  the  same  as  for  the  single-  channel 
case. 


4.  CONCLUSIONS. 

In  this  paper  we  have  developed  the  construction  of  optimum  receiver 
for  the  slowly  fluctuating  point  target  in  the  presence  of  white  or  colored 
noise.  It  is  concluded  that,  in  general,  in  the  case  of  the  n-channels  model, 
n2  matched  filters  (or  correlations)  are  required  to  obtain  the  receiver.  The 
matching  is  done  with  the  orthogonalized  signals  instead  of  the  transmitted 
signals.  The  optimum  receiver  computes  the  weighted  sum  of  the  outputs  of 
matched  filter.  These  weights  are  such  that  the  risk  in  detection  is  minimized. 
We  do  not  evaluate  the  performance  measures  of  these  optimum  detectors. 
However,  it  follows  from  the  distributional  properties  of  the  summands  that 
these  measures  can  be  derived  in  terms  of  weighted  sum  of  chi-squares.  When 
the  number  of  channels,  n,  is  large  these  performance  measures  can  be  ap¬ 
proximated  by  normal.  In  general,  approaches  discussed  in  [l]  and  [2]  are 
applicable. 
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